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Final  Report 

1.0  Introduction  and  Summary 

This  report  represents  the  Final  Report  on  Contract  N00014-90-C-0217,  BAA  90-06, 
“Applicability  of  Superconducting  Interconnection  Technology  for  High  Speed  ICs 
and  Systems”.  Additionally,  this  volume  includes  the  previous  four  quarterly  reports 
issued  as  part  of  this  effort. 

The  extremely  low-loss  characteristics  of  HTSC  interconnect  renders  it  an  attractive 
candidate  for  high  density-high  performance  digital  circuitry.  The  purpose  of  this 
program  was  to  investigate  these  properties  in  detail,  and  to  determine  the 
characteristics  of  systems  that  can  best  e^loit  these  properties.  To  this  end  a  study 
was  undertaken  to  more  fully  characterize  the  limitations  of  conventional 
interconnect,  in  particular  in  the  context  of  source-terminated  CMOS  devices,  a 
technology  well-adapted  to  ciyogenic  applications.  Thermal  and  mechanical 
considerations  of  HTSC-MCM  technology  were  explored  as  well  as  part  of  this 
effort.  Applying  Rent’s  Rule  to  system  connection  requirements,  it  was  determined 
that  HTSC  interconnects  can  significantly  reduce  the  required  levels  of  wiring 
and/or  increase  the  packing  density  of  conventionally  partitioned  IC’s  on  an  MCM. 
This  advantage  is  amplified  if  IC’s  are  first  optimized  for  MCM  application  (eg. 
using  larger  numbers  of  area-I/O  pads  at  the  chip  levels).  Commercial  applications 
of  the  technology  were  likewise  explored,  and  it  is  believed  that  the  development  of 
conventional  cryogenic-MCM  optimization  technology  could  prove  to  be  a  useful 
first  step  in  the  development  of  an  appropriate  HTSC-MCM  infra-structure.  In 
conclusion,  it  is  felt  that  HTSC-MCM  technology  has  the  potential  to  become  a 
high-payoff  commercially  viable  future  technology  for  high  performance  digital 
packaging. 
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1.1.  Background 

1.1.1.  High  Performance  Electronic  Packaging  Using  Multi-Chip  Modules 
with  Superconducting  Interconnects 

While  the  focus  of  this  program  is  on  the  applicability  of  superconducting 
interconnection  technology  for  high  speed  systems,  a  major  aspect  of  the  effort  is  on 
understanding  the  limitations  of  conventional  interconnect  materials  in  advanced 
packaging  applications.  Multi-chip  module  (MCM)  packaging  of  digital  electronic 
systems  not  only  allows  for  reductions  in  size  and  weight  by  using  bare  integrated 
circuit  (IC)  die  (instead  of  individually  packaged  circuits)  mounted  with  high 
density,  but  allows  for  higher  system  operating  speeds  as  well.  This  higher 
performance  is  made  possible  because  the  distances  between  the  various  IC  chips 
are  reduced,  which  correspondingly  reduces  the  time  delay  (propagation  delay) 
experienced  by  signals  going  from  one  part  of  the  ^stem  to  another  (the  velocity  of 
such  signals  being  limited  by  the  speed  of  light).  Since  the  clock  period  (or  time 
required  to  do  an  elementary  operation)  must  exceed  the  sum  of  this  signal 
propagation  delay  plus  the  time  it  actually  takes  the  IC  chips  to  perform  their  logic 
functions,  the  system  performance  can  only  be  increased  if  both  the  logic  and 
interconnect  delays  are  reduced.  In  recent  years,  improvements  in  semiconductor 
technology  have  sharply  increased  IC  speeds,  but  it  has  become  increasingly  difficult 
to  achieve  the  full  measure  of  this  benefit  in  higher  system  speeds  because  of 
packaging  limitations. 

The  highest  density,  and  hence  shortest  interconnect  delays,  are  achieved  when  an 
MCM  substrate,  large  enough  to  hold  all  of  the  chips  required,  is  essentially  "tiled" 
with  IC  die,  mounted  as  closely  together  as  the  bonding  technology  allows.  This  type 
of  large,  very  high  density,  "system  on  an  MCM"  module  is  illustrated  in  Figure  1.1. 
with  such  maximal  density  MCM  packaging,  a  single  6"  MCM  substrate  might  hold 
500  or  more  VLSI  chips,  allowing  a  complete  very  high  performance  system  to  be 
placed  on  a  single  MCM.  Unfortunately,  it  is  extremely  difficult  to  fabricate  such  an 


Applicability  of  Superconducting  Interconnection 
Te^ology  for  Hi^  Speed  ICs  and  Systems 


DARPA  90-06 


1-2 

October  1992 


^  o 

-3  « 

Q  C  O) 

.00 
c\J  .9  C 
E 

•  -  o  *0 
w  (D 

C  ^  CD 

0  Q. 

Q  o  CO 
-c  0  ^ 
.95  fc  .g> 
X  <  X 


FINAL  REPORT 


MCM  using  conventional  metal  (eg.,  copper)  interconnects  because  the  resistivity  of 
normal  metals  are  too  high.  For  this  6"  MCM  example,  over  100,000  signal 
interconnect  wires  would  be  needed  to  make  the  connections  between  the  VLSI 
chips,  requiring  a  total  wire  length  of  nearly  two  miles  on  the  MCM!  Actually,  the 
lithography  used  in  IC  fabrication  is  capable  of  producing  this  much  wire  with  only 
two  layers  of  1.5/im  wide  signal  lines  (on  a  6jum  pitch),  but  with  normal  metal  lines, 
their  resistance  would  be  hopelessly  too  high  (4600S2  for  a  12"  copper  line,  whereas 
it  should  be  50Q  or  less  for  proper  signal  propagation).  With  normal  metals,  the 
only  way  to  keep  the  line  resistance  low  enough  is  to  make  the  wires  bigger  in  area, 
which  reduces  the  density  of  the  wires  in  each  signal  interconnect  layer, 
necessitating  the  use  of  very  large  (unpractically  large  from  a  fabrication  standpoint) 
number  of  metal  interconnect  layers.  For  example,  for  the  6"  MCM  case  discussed, 
over  60  copper  signal  wire  layers  would  be  required. 

High  temperature  superconductors  (HTSCs)  offer  an  exciting  solution  to  this  high 
density  MCM  dilemma.  These  materials,  when  cooled  to  cryogenic  temperatures  of 
the  order  of  TT’K,  have  essentially  zero  resistivity,  and  hence  could  in  principle  be 
used  to  wire  such  maximal  density  MCMs  with  only  two  signal  interconnect  layers. 
Not  only  would  this  make  fabrication  of  such  large  "system  on  a  substrate"  MCMs 
with  their  veiy  short  interconnect  signal  propagation  delays  practical,  but  there 
would  be  a  matching  reduction  of  the  logic  delay  as  well.  Because  the  velocities  and 
mobilities  of  the  electrons  or  holes  in  the  semiconductor  devices  from  which  the  ICs 
are  built  are  greatly  increase  at  low  temperatures,  typically  the  logic  delays  are 
reduced  by  a  factor  of  about  2.S  when  the  chips  are  operated  at  the  same 
temperatures  as  would  be  used  for  HTSC  interconnections. 

Implementation  of  a  digital  electronic  system  in  an  HTSC  MCM  packaging 
approach,  through  the  balanced  reduction  of  logic  delay  (from  low  temperature  IC 
operation)  and  signal  interconnect  propagation  delay  (from  maximal  density 
packaging),  offers  a  major  increase  in  system  performance  (eg.,  2.5  time  faster)  over 
normal  systems  using  the  same  type  of  ICs.  In  other  words,  if  a  computer  could  run 
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at  an  lOOMHz  clock  rate  normally,  then  250MHz  should  be  achievable  using  HTSC 
MCM  packaging,  which  would  greatly  increase  its  value  (eg.,  if  it  were  worth  $1M 
with  conventional  packaging,  it  would  be  worth  about  $2.5M  in  the  HTSC  MCM 
version).  Of  course,  in  many  military  applications  the  fact  that  the  size  and  wight 
have  been  dramatically  reduced  (to  a  single  6"  x  6"  substrate)  is  of  crucial 
importance,  over  and  above  the  fact  that  one  faster  processor  can  do  the  work  of  2 
to  4  normal  ones  (which  further  save  size,  weight,  cost  and  power). 

While  the  advantages  of  the  use  of  HTSC  intercoimects  in  digital  MCMs  are  obvious 
and  impelling,  substantial  research  and  technology  development  effort  will  be 
required  to  achieve  this  result.  The  signal  current  densities  in  the  fine  (eg.,  W  = 
1.5;im  or  ^m )  superconducting  signal  lines  are  achievable  at  present  only  with  high 
quality  thin  film  HTSC  materials.  These  are  generally  grown  on  crystalline 
substrates  in  a  single  flat  layer,  whereas  the  MCM  application  requires  the  growth  of 
several  patterned  layers  with  thick  layers  of  reasonably  low  dielectric  constant 
insulators  (crystalline  or  amorphous)  in  between  them,  as  illustrated  in  Figure  12. 
In  addition  to  these  basic  HTSC  materials  growth  issues,  it  will  be  necessary  to 
develop  a  practical  HTSC  MCM  fabrication  process  for  making  HTSC  MCM 
structures  like  that  shown  in  Figure  12,  including  not  only  growing  and  patterning 
the  HTSC  signal  and  power  and  ground  planes  separated  by  the  appropriate 
dielectric  layers,  but  handling  the  interlayer  via  connections,  die  attach  and  lead 
bonding,  supply  bypassing,  I/O  connection,  etc.  problems  as  well.  In  addition,  there 
needs  to  be  associated  with  the  technology  development  effort  a  strong  electronic 
systems  presence  to  provide  focus  to  insure  that  the  HTSC  MCM  technology 
developed  will  be  appropriate  for  realistic  digital  systems  applications. 

The  needs  of  this  HTSC  MCM  technology  development  effort  extend  from  basic 
HTSC  materials  research  through  fabrication  process  development  to  the  creation 
of  a  business  base  capable  of  commercially  supplying  such  a  technology  for  a  wide 
variety  of  high  performance  digital  systems  applications,  both  military  and 
commercial.  This  will,  of  necessity,  require  a  team  effort  to  accomplish  these  ends. 
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with  involvement  of  imiversities  (for  basic  HTSC  materials  research),  merchant 
HTSC  companies  (within  which  to  develop  the  HTSC  MCM  technology  to  insure  it 
becomes  commercially  available  to  all  US  companies),  and  companies  capable  of 
furnishing  required  supporting  technologies  (eg.,  cryogenic  coolers  and  packaging, 
ICs,  testing,  CAE  tools,  etc.),  plus  an  electronic  systems  house  to  maintain  the 
proper  applications  focus  to  the  effort. 

The  payoff  from  such  an  HTSC  MCM  technology  development  effort  should  be 
enormous.  Not  only  would  such  commercial  digital  electronics  systems  as 
computers  (eg.,  supercomputers,  mainframes,  mini-supercomputers  and 
super-minicomputers),  communications  systems  and  specialized  signal  processors 
such  as  pattern  recognition  or  medical  imaging  systems  be  potential  applications, 
but  a  wide  variety  of  military  signal  processing  (eg.,  RADAR,  ELINX  ECM,  "smart 
weapons",  etc.),  communications  and  other  military  system  applications  could 
greatly  benefit  fi-om  this  exciting  technology  as  well.  In  addition  to  MCM  substrate 
applications,  the  same  HTSC  technology  to  be  developed  for  HTSC  MCMs  could  be 
applied  to  many  other  applications  as  well.  For  example,  this  HTSC  MCM 
technology  could  be  utilized  to  implement  an  HTSC  backplane  which,  in  only  a  6" 
width,  could  achieve  a  2.5  Terabits/second  raw  data  transfer  rate  at  clock  rate  of  only 
lOOMHz,  or  up  to  250  Terabits/second  if  the  ICs  were  fast  enough  to  handle  a  lOGHz 
clock  rate  (that  rate  would  be  no  problem  for  the  fineline  HTSC  interconnects). 
This  is  only  one  of  many  other  possible  applications  for  this  exciting  technology. 
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1.2.  Commeclal  Applications 

We  began  this  program  knowing  that  if  HTSC-MCM  technology  coiild  be  developed 
and  the  benefits  of  lower  temperature  operations  are  realized,  the  MCM  industry 
would  still  encounter  difficulties  having  this  technology  accepted  by  the  high  volume 
commercial  community.  This  part  of  the  study  was  aimed  at  understanding  the  high 
volume  market,  both  commercial  and  military,  and  determine  the  criteria  under 
which  HTSC-MCM  would  be  accepted  by  system  producers. 

About  a  year  ago,  having  met  with  a  number  of  commercial  systems  houses,  we 
realized  that  while  HTSC-MCM  technology  held  great  promise,  much  development 
is  still  needed  to  achieve  the  ultimate  goals  of  the  program.  However,  as  we  got  into 
it  in  more  depth  we  discovered  that  the  HTSC-MCM  fovmdries  could  become 
unstable  since  their  financial  backers  may  not  have  either  the  patience  and  the 
stamina  to  wait  for  this  technology  and  its  market  to  develop.  We  subsequently 
recommended  that  this  program  should  be  modified  to  allow  those  foundries  to 
consider  the  possibility  of  fabricating  conventional  MCMs  but  operating  them  at 
ciyotemperature.  This  approach  has  several  benefits:  (a)  Cold  CMOS  does  exhibit 
significant  performance  improvement  below  -100°C;  (b)  Some  companies  are 
already  considering  the  possibility  of  cooling  their  systems  to  gain  performance, 
meaning  a  near  term  market  could  be  addressed  by  the  HTSC  foundries;  (c)  These 
foundries  need  to  develop  the  cryotechnology  anyway,  including  cryogenic 
assemblies,  cryotesting  and  packaging,  and  IC  optimization  for  cryotemperature 
operation.  We  consequently  crusaded  the  industry  infrastructure  to  alter  the 
HTSC-MCM  program  and  make  it  slightly  more  general  The  new  program  that  we 
recommended  should  really  be  cryoelectronics  which  includes  HTSC-MCM,  IC 
optimization  and  cryogenic  cooling,  etc.  This  program  received  the  support  of 
industry,  universities,  and  DARPA.  This  program  indeed  becomes  more  important 
in  light  of  the  finding  by  SIA  that  within  ten  to  fifteen  years  the  IC  industry  will  face 
major  challenges  if  they  plan  to  continue  the  same  rate  of  on-chip  performance  and 
complexity.  Not  only  is  the  cost  of  new  CMOS  foundries  reaching  exhorbitant  levels 
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($1B  per  foundry)  but  now  they  have  to  rely  on  new  processing,  testing  and 
packaging  technologies  to  reach  their  goals.  The  result  is  an  unquestionable  high 
cost  in  R&D  and  a  significant  delay  to  the  marketplace.  Cryoelectronics  could 
relieve  a  great  deal  of  the  pressure. 

More  recently  however  we  reached  yet  another  plateau.  It  became  clear  to  us  that 
even  conventional  MCM’s  are  having  difficulty  being  accepted  by  system  houses. 
We  believe  the  reason  for  this  to  be  that  MCM’s  are  being  populated  bv  IC’s 
optimized  for  single  packages  and  not  for  MCM’s.  Initial  studies  RC  Eden  and 
Doug  Paffel  clearly  showed  that  if  IC’s  are  optimized  for  MCM’s  new  and  major 
strides  can  be  achieved  in  the  near  term.  The  following  is  a  summary  of  our  findings: 

The  MCM  technology  brings  in  two  major  features:  (1)  Wiring  densities  comparable 
to  those  used  on  chips,  and  (2)  Via  pitch  transformer  capability  that  can  couple  the 
IC  to  the  PCB.  With  these  two  features  we  now  can  do  the  following:  (a)  We  can 
wire  chips  to  each  other  using  extremely  high  numbers  of  I/O  contacts.  This  means 
that  in  going  from  one  chip  to  the  other,  the  connectivity  does  not  go  through  major 
discontinuity.  Interchip  intercoimect  will  appear  much  like  the  on  chip  wiring. 
Therefore  if  the  IC’s  have  area  pad  I/O  we  could  obviate  the  dis-economic  urge  to 
design  very  large  IC’s.  Since  large  IC’s  are  designed  to  circumvent  the  constraints  of 
low  peripheral  IC  I/O  counts,  the  combination  of  IC  area  I/O  and  high  wiring  density 
MCMs  could  give  the  electronics  industry  the  ability  to  wire  many  small  area-I/O 
chips  on  MCMs  and  make  them  look  as  if  they  are  a  single  jumbo  chip.  The  benefits 
of  this  approach  are  indeed  staggering:  (a)  By  using  small  chips  we  continue  to 
operate  on  the  high  yield  part  of  the  yield  curve,  this  means  cost  per  gate  can  be 
reduced  significantly,  (b)  Functionality  that  was  compromised  out  to  the  peripheral 
I/O  limitation  can  be  reinstated  adding  an  additional  factor  in  performr  p..  (c) 
Using  smaller  blocks  of  chips  gives  us  higher  flexibility  of  using  these  chips  thereby 
adding  to  the  utility  and  perhaps  product  life  cycle  and  (d)  Earlier  market 
penetration  since  smaller  chips  can  then  be  produced  more  quickly. 
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Having  said  all  of  the  above  it  now  becomes  clear  that  devices  optimized  not  only  for 
ciyotemperature  operation  but  also  for  MCM  packaging  could  lead  to  performance 
improvements  exceeding  16x  and  price/performance  ratios  possibly  exceeding  64x. 
It  is  for  this  reason  that  now  we  recommend  that  the  HTSC-MCM  program  should 
be  broadened  to  include  MCM  and  cryotemperature  optimization  ranging  from 
room  temperature  to  80°k.  This  program  should  give  system  houses  great  incentives 
to  use  MCMs  over  a  wide  range  of  temperature  operation. 
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1.3.  Source-Terminated  Lossy  Line  Transient  Response 
Characteristics 

The  principal  integrated  circuit  technology  of  interest  for  cyoelectronic  applications 
(with  or  without  HTSC  interconnects)  is  silicon  CMOS.  Far  and  away  the  most 
common  signal  interconnection  scheme  employed  with  CMOS  ICs  is  the  source 
terminated,  open-load  transmission  line  configuration  shown  in  Figure  1.3.  While 
the  Fourier  transform  approach  (see  Section  2  and  3)  for  very  accurately  modeling 
the  transient  response  of  small  lossless  dielectric  lossy  conductor  transmission  line 
interconnects  (eg.  as  used  in  MCMs)  will  work  for  any  signal  interconnect  system 
having  reasonably  linear  driver  and  load  characteristics,  our  focus  thus  far  has  been 
principally  on  the  CMOS-like  source  terminated  open  load  configuration  of  Figure  13. 

Before  presenting  the  results  of  the  detailed  computer  modeling  of  these 
interconnects,  it  is  worth  spending  a  few  minutes  effort  in  trying  to  understand  from 
an  analytical  viewpoint.  The  lossless  =  0  and  -  0)  case  for  the  source 
terminated  open  load  (STOL)  configuration  of  Figure  13  is  discussed  in  the  last 
quarterly  report  (Section  2  of  this  document),  and  should  best  be  read  before 
proceeding  with  this  final  report  material.  Note  that  the  material  in  the  last  three 
quarterly  reports  (Sections  3.2,  3.3  and  2  of  this  document)  is  presented  as  one 
continuous  text,  with  figure  numbers,  equation  numbers  and  reference  numbers 
consistent  between  them  (only  the  reference  number  continuity  is  maintained  in  this 
final  report  section). 

As  shown  in  Figures  36  through  49  in  Section  2,  for  the  lossless  (R^  and  R^  =  0) 
conductor  case,  the  lossless  STOL  load  end  voltage  response,  ^2,  is  a  sequence  of 
output  steps,  beginning  at  1  ^  (where  is  the  one-way  preparation  delay)  each 

step  of  duration  2  (the  round  trip  propagation  delay).  Equations  49  and  51  give 
the  amplitudes  of  the  first  two  steps  (^^  <  t  <  3t^  and  3^  ^ 

respectively).  The  major  pulse  distortion  effect  which  is  introduced  by  lossy 
conductors  is  caused  by  the  dc  series  resistance  (primarily  of  the  center  conductor  in 
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coax,  as  the  larger  ground  plane  is  lower  resistance),  in  Figure  1.3.  The  "zero 
order"  approximation  to  the  lossy  line  response  is  obtained  by  noting  in  Figure  1  that 
at  t  there  is  a  current  flowing  through  a  series  resistance  distributed 
along  the  line,  and  so  the  output  voltage  will  be  reduced  by  an  amount  AF  =  R^ 
below  that  given  in  Equation  49.  This  means  that  to  "zero  order",  the  initial  load 
voltage,  Fj,  should  be  given  by: 


=  V) 


2F  Z 

in  o 


^in 


+  z. 


1 


Eq.1.1 


Immediately  after  ^  a  negative  reflected  current  pulse  in  Figure  1.3)  starts 
back  toward  the  source,  which  progressively  reduces  the  (/^  -  ohmic 

conductor  drop  on  the  line.  To  "zero  order",  this  conductor  drop  goes  to  zero  when 
the  pulse  gets  back  to  the  source,  which  is  seen  at  t  =  2 at  the  source,  but  not 
until  t  -  3  at  the  load  end.  Hence,  the  "zero  order"  at  the  load  end  we  see  at 
/  =  3  (the meaning  "just  before  the  reflection  discontinuity  arrives") 


^in  + 


Eq.l2 


which  is  the  same  as  the  lossless  step  height  of  Equation  49.  In  between  <  =  and 
t  =  3  ^  the  ouft)ut  voltage  will,  to  "zero  order"  just  go  linearly  between  the  values 
given  by  Equations  1.1  and  1.2,  since  the  I2  pulse  moves  at  constant  velocity, 
"erasing"  the  (J^  -  ohmic  voltage  drop  linearly  in  time. 
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These  "zero  order"  expressions  are  reasonably  accurate  when/?^  is  small  but  do  not 
count  the  fact  that  the  current  along  the  line  must  decrease  somewhat  as  the  IR^ 
drop  reduces  the  voltage  (since  the  ratio  of  voltage  to  current  is  Z^).  This  simply 
implies,  of  course,  an  exponential  decay  of  the  initial  pulse  height  with  /2Z^ 
instead  of  the  linear  decay  of  Equation  1.1.  Hence,  to  first  order  (as  shown  for 
R.^  =  Z^in  Equation  7,  Section  3.2)  we  have  the  output  voltage  at  r  =  given  for 
(only)  loss  by: 


= 


)  = 


2Z. 


for  R^^  (only)  Eq.  13 


The  various  current  reflections  which  make  Equation  13  different  fi’om  eq.  1.1  do 
not  all  die  out  completely  at  r  =  3  so  that  (simply  updating  Equation  8  for  the 
R.^  ^  general  case)  we  will  have  at  r  =  3  t“, 


^2  = 


2V.  Z 

in 

Rin  +  Z„ 
in  o 


-0.48 


2Z. 


2Z. 


for  7?^  (only)  Eq.  1.4 


In  the  lossless  and  i?^^(only)  lossy  transient  response  characteristics,  the  step 
response  is  marked  by  discontinuities  in  output  voltage,  V^zX  t  =  » 3  ^ » 5  ^ , 

etc.  These  "zero  risetime"  output  steps  are  an  artifact  of  the  zero  risetime  of  the 
input  and  the  infinite  bandwidth  assumed  for  the  transmission  line  because  we  have 
ignored  the  skin  effect  losses  on  the  conductors  (primarily  the  center  conductor  in 
coaxial  lines).  When  skin  effect  loss  is  included  (eg.,  see  Subsection  1.23,  Pages  1-9 
to  1-13  and  Figure  1  of  Section  33,  Equations  9-23)  there  abrupt  transitions  are 
shown  to  have  a  finite  risetime  and  a  complementary  error  function  shape 
(Equations  14  and  15  and  Figures  2  through  6).  A  practical  approach  for  including 
both  the  Rj^  and  skin  effect  loss  characteristics  in  an  analytical  model  was  given  in 
Subsection  1.2.5,  Pages  1-15  through  1-17  and  Figures  7  through  15. 
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These  analytic  approximations  handle  reasonably  low  loss  lines  fairly  accurately 
except  in  the  "tail"  region  (approaching  asymptotic  response),  where  small  errors 
tend  to  accrue.  This  is  also  the  case  for  the  Quad  Design  XNS  tools  (which  is  quite 
satisfactory  for  the  engineering  design  application  for  which  these  tools  are 
intended,  where  a  1-2  percent  error  in  the  far  tail  response  is  not  of  consequence  in 
digital  circuits).  However,  for  highly  lossy  and  very  small  lines,  some  of  the 
approximations  in  standard  analyses  are  inadequate  and  place  the  accuracy  of  the 
resulting  transient  responses  in  question.  The  intent  of  the  effort  described  in  detail 
in  Section  3.3  was  to  model  these  extreme  cases  with  a  high  degree  of  accuracy  (but 
limited  to  linear  drivers  and  not  having  as  fast  run-times  as  needed  for  commercial 
engineering  tools).  Because  skin  effect  can  be  easily  modeled  in  the  frequency 
domain,  the  analysis  is  done  initially  in  the  frequency  domain,  and  then  an  inverse 
Fourier  transform  (8192  point)  is  carried  out  to  produce  the  time  domain  impulse 
response,  which  is  subsequently  integrated  to  give  the  step  response.  While  small 
copper  coax  lines  were  modeled  in  this  work  so  that  analytical  expressions  (using  the 
her  and  bei  Kelvin  functions  given  in  Equations  40  and  41)  for  the  ac  series  resistance 
and  internal  inductance  of  the  center  conductor  are  known,  other  geometries  could 
be  modeled  using  the  approximations  given  in  Section  2.4  of  Reference  5.  In  the 
analysis,  both  the  characteristic  impedance,  Z^,  and  the  propagation  constant,  y,  are 
taken  as  complex,  frequency  dependent  quantities.  The  analysis  was  done  using  the 
Bimillenium  "HiQ"  engineering/scientific  software  package,  running  wither  on  a 
Macintosh  fx  or  a  Powerbook  170  computer.  The  step  response  results  are  stored  as 
an  8192  point  vector  and  are  analyzed  and  plotted  using  Microsoft  EXCEL. 
Because  EXCEL  4.0  cannot  handle  more  than  4000  points  in  a  plot,  it  was  necessary, 
in  order  to  plot  "full  time  span"  results  like  those  in  Figures  1.4, 1.6, 1.9  etc.,  to  create 
a  second  4096  point  output  file  (skipping  every  other  data  point  in  the  8192  point 
file).  The  expanded  time  scale  plots  (Figures  1.5, 1.7, 1.10,  etc.)  were  made  from  the 
original  8192  point  files  for  the  best  time  resolution  (4ps  in  Figures  1.4  through  1.16 
and  2ps  in  Figures  1.17  through  1.26,  except  for  Figmes  1.6  through  1.8  which  are  a 
2.Sps  step  size). 
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1 .4.  Small  Lossy  Coax  Transient  Response  Simulation  Results 
for  Fixed  Overshoots 


A  key  question  in  MCM  technology  is  what  is  the  maximum  allowable  line  resistance 
(dc  or  ac)  for  the  signal  interconnects.  For  example,  in  References  5  and  6  a 
^MAX  value  was  taken  for  the  maximum  conductor  resistance  for  most  cases, 
or  =  25Q  considered  for  so  called  "self-terminated  lines"  (see  discussion  by 
Equations  SO  and  51).  One  question  is  up  to  what  frequency  the  conductor 
resistance  should  be  kept  to  signals  of  a  given  risetime.  Another  question 

is  how  effectively  the  distributed  ohmic  dc  or  skin  effect  resistance  can  be 
compensated  by  reduction  in  (the  source  termination  resistance  at  the  driver). 
Also,  what  would  happen  to  the  voltage  waveforms  if  all  output  drivers  had  the  same 
reduced  R^  values,  even  if  they  were  driving  shorter,  lower  resistance  lines.  The 
answers  to  these  questions  allow  us  to  determine  how  small  the  copper  signal  lines 
may  be  made  for  a  given  size  and  performance  level  MCM. 


In  Section  2, 12"  long  300“K  copper  coax  lines  (eg.,  longest  lines  on  a  6"  MCM)  were 
analyzed  with  center  conductor  diameters  of  1.0  mils,  0.5  mils  and  0.25  mils,  for  a 
range  of  R.^  values  (eg.,  30Q ,  35Q,  40Q,  45Q,  50Q  and  55S2  in  Figures  50  and  51). 
Also,  one  6"  line  length  case  with  a  0.5  mils  center  conductor  diameter  and  a  similar 
range  of  values  was  presented.  In  this  section,  simulation  results  are  presented 
for  a  full  matrix  of  line  lengths  and  center  conductor  diameters.  Specifically, 
diameters  of  1.0  mils,  1/VJ  mils  (0.7071  mils),  0.5  mils,  l/2>^  mils  (03535  mils)  and 
035  mils  are  analyzed  in  line  lengths  of  12.0"  (30.48  cm)  and  6.0"  (15.24  cm).  The 
V?  stepping  of  center  conductor  diametei^  was  done  to  give  a  factor  of  2  stepping  in 
R^  between  the  various  cases. 

The  other  major  difference  in  the  simulations  of  Figures  1.4  through  1.26  is  the  way 
in  which  the  driver  source  resistance,  R.^  was  selected  for  the  various  curves.  The 
focus  is  on  illuminating  how  helpful  the  self-termination  approaches  can  be  in 
allowing  (by  reducing  R^ )  lines  that  would  otherwise  have  unacceptably  poor 
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response  (due  to  excessive  series  resistance)  to  be  used  in  digital  applications.  In 
general,  good  digital  design  has  some  specified  "budget"  for  line  aberrations  such  as 
overshoot,  undershoot,  "ringing",  etc.,  just  as  it  has  for  crosstalk  and  other  factors 
which  eat  into  the  total  noise  margin  budget.  A  typical  number  might  be  a  5%  or 
10%  allowance  for  overshoot,  depending  on  the  IC  technology.  Hence,  in  these 
figures,  curves  are  presented  for  selected  R.^  values  which  give  peak  overshoot 
values  of  +20%,  + 10%,  and  +5%  (as  defined  in  Equation  50),  along  with  one  or 
more  curves  at  near  zero  overshoot.  The  latter  is  difficult,  because  in  the  <  1% 
overshoot  region,  a  lot  of  subtitles  in  the  accurate  transient  waveforms  arise 
(principally  due  to  the  complex,  frequency  dependent  )  which  make,  for  example, 
the  R.  for  an  overshoot  of  +  0.6%  quite  different  than  the  R.  to  give  a  0.05% 
overshoot.  (This  is  because  the  point  of  peak  overshoot  moves  away  from  ^  =  3  ^ 
as  assumed  for  the  lossless  case  in  Figures  37  and  38).  You  can  see  these  effects  in 
Figures  1.11, 1.14  and  1.25  particularly. 

For  the  lossless  case  (eg.  Equations  49  and  50  or  Figures  37  and  38),  for  a  nominal 
Z  =  50112,  an  overshoot  of  0%  would  be  obtained  with  R.  -  50Q,  +5%  with  R.  = 
45.238095238Q,  + 10%  overshoot  with  i?.  =  40.90909090  and  +20%  with  = 
33.3333330.  For  cases  with  low  values  of  R^,  such  as  the  6.0"  1.0  mil  center 
conductor  diameter  case  of  Figures  1.17  and  1.18  {R^  =  6.0980),  R^  = 
34.161770  gave  +20%  overshoot,  R.  =  42.655440  gave  +10%,  and  R.  = 
47.521960  gave  +5%  overshoot.  These  i?^values  are  close  to  those  predicted  for 
the  lossless  case.  This  is  not  the  case  for  near  zero  overshoot  For  example,  a  value 
of  R.  -  53.0089240  gives  an  overshoot  of  +0.962%  (at  x  =  3.186ns  and  R.  = 
62.0000  gives  an  overshoot  of  +  0.01315%  (at  x  =  7.450  ns),  whereas  in  the  lossless 
analysis  not  including  the  complex  Z^  any  R^  over  500  would  be  underdamped 
(negative  overshoot). 

IWo  different  sets  of  analysis  runs  were  made  for  the  12"  long,  1.0  mil  center 
conductor  diameter  case,  the  first  (Figures  1.4  and  1.5)  at  a  ~  4ps  step  size 
(after  ifft),  and  the  second  (Figures  1.6  and  1.7,  summarized  in  Figure  1.8)  at  a  2.5ps 
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Step  size.  The  results  were  virtually  identical  for  the  two  analyses.  For  example, 
with  =  41.7886594Q,  the  2.5ps  run  gave  an  overshoot  of  10.0000817%  at  6.040ns 
and  the  4ps  run  gave  9.9987777%  at  6.040ns,  only  a  20.4//y  difference  in  peak  output 
in  a  1  volt  pulse  amplitude.  Additional  verification  of  the  accuracy  of  the  calculations 
is  the  fact  that  there  is  virtually  zero  output  seen  at  prior  to  the  nominal  pulse 
arrival  time,  (where  =  c/Ve^  =  1.518  x  10®meters/s  for  the  =  3.90 

dielectric  constant  assumed,  or  1/v^  =  1673ps/inch).  For  /  =  12",  =  2.007ns,  or 

for  a  /  =  6"  line  length,  =  1.00ns.  Note  in  the  most  expanded  time  scale  /  =  12" 
plot.  Figure  1.7,  the  extrapolated  pulse  arrival  time  is  indeed  2.0075ns  (the  2.0075ns 
step  showing  zero  output  and  the  next  time  point  at  2.010ns  showing  finite  output). 
Another  measure  of  the  accuracy  of  the  calculation  is  the  fact  that  while  the  original 
calculation  is  done  in  the  frequency  domain,  (inverse)  Fourier  transformed  to  the 
time  domain  impulse  response,  with  the  result  integrated  to  give  the  step  response, 
without  any  amplitude  or  slope  correction  the  asymptotic  pulse  amplitude  goes 
precisely  to  1.00000000  volt.  This  accuracy  is  due  in  part  to  the  numerical  precision 
of  the  Bimillenium  HiQ  software,  coupled  with  care  in  getting  the  physics  into  the 
problem  correctly  (ie.,  being  very  careful  that  the  frequency  domain  pin-phase 
characteristics  do  not  violate  the  Kramers-Kronig  relationship)  and  avoiding 
approximations  which  would  induce  artifacts  into  the  output. 

As  mentioned  previously,  in  digital  systems  it  is  a  common  design  approach  to  assign 
a  specific  overshoot/undershoot  budget  as  part  of  the  overall  noise  margin  budget. 
This  differs  from  the  peak  overshoot  value  in  that  it  represents  the  actual  signal 
error  at  the  time  the  logic  decision  is  made.  The  peak  overshoot  value  may  lead  to 
improper  circuit  function  or  even  degraded  IC  reliability.  The  error  budget 
concerns  the  accuracy  of  the  decision  process  (in  analogy  to  the  signal  settling  time 
in  an  analog  circuit  or  A/D).  As  seen  in  Figures  1.4  through  1.26,  in  this 
source-terminated,  open-load  configuration  (Figure  3),  substantial  improvements  in 
signal  risetime  at  the  end  of  lossy  lines  can  be  achieved  by  reducing  the  driver 
resistance,  R.  below  the  normal  R.  =  Z  value  to  compensate  for  the  ohmic  line 
resistance.  It  is  important  to  note,  however,  that  reducing  R^  may  in  fact  lengthen 
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the  "settling  time"  necessary  to  get  the  signal  amplitude  within  a  specified  error 
range.  Some  of  these  timing  relationships  are  illustrated  for  the  12"  1.0  mil  center 
conductor  case  in  Figure  1.8,  and  for  the  12"  0.5  mil  case  in  Figure  1.13.  The  most 
common  risetime  specifications  are  for  10%  to  90%  or  20%  to  80%  of  final  signal 
amplitude  values.  For  the  12"  1.0  mil  case,  the  10%  to  90%  risetime  value  with  0% 
overshoot  is  1387ns,  or  69%  of  the  nominal  delay  time,  which  value  is  halved  to 
0.6588ns  (32.8%  of  ^ )  by  reducing  R.^  to  46.9Q  to  give  a  5%  peak  overshoot,  or 
even  cut  to  205.4ps  if  is  reduced  to  32.8Q  for  a  20%  overshoot.  The  problems 
with  this  latter,  20%  overshoot  case,  noted  in  the  row  with  the  [110%  -  90%]  last 
crossing  amplitude  entry  in  the  left  column,  even  if  the  decision  error  budget  were 
±10%,  the  "setting  time"  would  be  6.3125ns,  and  because  of  the  ringing  undershoot, 
if  the  budget  were  a  tighter  2%  value  (shown  in  the  [102%  -  98%]  row  for  the  +  20% 
column),  we  would  have  to  wait  10.62ns  for  the  signal  to  recover  to  above  98%  of 
final  value  (as  can  be  seen  in  Figure  1.4).  Such  tight  error  budgets  (eg.  1%  to  2%  or 
less)  are  ^ically  associated  with  A/D  conversion  as  threshold  decision  circuits, 
values  of  5%  or  possibly  10%  would  be  more  common  in  all  digital  designs. 
However,  remember  that  this  12"  1.0  mil  case  with  a  total  dc  line  resistance  of 
122Q  would  be  considered  a  good  high  quality  signal  line  in  MCM  use.  Any  shorter 
lines  which  used  the  same  R.^  value,  as  noted  previously,  would  have  larger  peak 
error  budgets. 

Of  particular  interest  is  the  case  of  lines  in  the  "self-terminating"  region,  with  R^ 
values  not  small  in  comparison  to  Z^.  Figure  1.13  shows  this  same  detailed  timing 
data  for  the  case  of  a  12"  long  line  with  a  03  mil  center  conductor  diameter  which 
has  R^^  =  48.78Q.  As  seen  in  Figure  52,  with  a  nominal  R.  -  Z  =  50Q  value,  the 
10%  to  90%  risetime  is  a  dismal  5.3ns  (2.6  times)  the  propagation  delay.  If  we 
reduce  to  42.7Q  for  a  0%  overshoot,  the  10%  -  90%  risetime  improves  to  2.417ns 
(still  1.2  times  t^^,)  or  if  we  further  reduce  R.^  to  36.12S2  for  a  5%  peak  overshoot  we 
can  get  a  1.68ns  risetime  (still  84%  of  ^^).  Note  that,  as  seen  in  Figure  38,  shorter 
lines  with  this  R.^  could  have  up  to  a  16%  overshoot,  which  could  well  be  acceptable. 
If  more  radical  reductions  in  R.  are  made,  such  as  to  R.  =  303Q  for  a  10% 

i/I  I/I 
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overshoot,  the  10%-90%  risetime  (1.07m)  approaches  half  of  ^  but  the  peak 
overshoot  on  shorter  lines  with  this  would  approach  25%  which  might  not  be 
acceptable.  As  seen  in  Figure  1.12,  the  risetime  can  be  reduced  to  487ps  by 
dropping  R.^  to  19.842Q,  giving  a  20%  peak  overshoot,  but  if  this  R.^  value  were 
used  for  shorter  lines,  the  peak  overshoot  could  be  severe,  up  to  43%  (Equation  12) 
with  heavy  ringing  (eg.,  like  Figure  48). 

While  the  risetime  degradation  is  serious  and  pulse  overshoot  can  be  severe  if 
comtant  R.^  is  used  on  all  drivers,  clearly  (Figures  1.9  through  1.13)  reduced  R^, 
"self  terminated  line"  interconnect  approaches  could  be  useable  for  some  digital 
applicatiom  in  the  R^^  =  25Q  to  50Q  range  (R^  =  ZJl  to  Z^),  using  R^  values  in 
the  25Q  to  40Q  range.  If  the  wire  diameter  is  reduced  another  factor  of  2  however, 
as  shown  in  the  12"  0.25  mil  center  conductor  diameter  case  of  Figure  1.16,  even 
reducing  R.^  to  zero  will  not  bring  this  =  195Q  case  to  near  critical  damping. 
Even  with  an  R.^  =  30Q  source  resistance,  the  90%  point  of  signal  respome  in  not 
reached  until  nearly  8m  after  nominal  pulse  arrival  (4  times  ^  ).  As  seen  for  /  = 
6.0"  in  Figures  125  and  1.26,  in  shorter  lines,  such  small  (1/4  z^)  copper  conductor 
diameter  offers  some  hope,  but  even  a  =  25Q  value  will  not  bring  the  r  =  3 
value  up  to  a  positive  overshoot,  and  a  negative  drive  resistance  (R^  =  -IJIOIQ, 
which  would  presumably  oscillate  short  lines  would  be  required  to  reach  a  +20% 
overshoot  for  this  R^^  =  97.6Q  case.  However,  for  short  interconnects,  even 
underdamped  respome  characteristics  may  be  useable  if  IC  speeds  are  not  very  fast, 
as  the  absolute  magnitudes  of  the  interconnect  delays  are  reduced.  Note,  however, 
that  we  have  ignored  the  effects  of  load  capacitance,  C^,  here,  and  the  simple  {R^  + 
Rd^Ct  time  comtant  effects  can  be  severe  for  small  normal  metal  lines. 

In  summary,  these  calculatiom  of  small  normal  metal  lines  respome  characteristics 
show  that  the  assumptiom  for  Rj^j^x  HTSC  MCM/normal  metal 

interconnect  MCM  tradeoff  studies  (References  5  and  6)  were  valid  (jmt  as  the 
detailed  "real  world"  netlist/wireability  study  of  References  9  showed  that  our 
statistical  wiring  model  assumptions  were  valid).  These  studies  show  that  as  the 
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^  number  of  die  on  MCMs  increases,  and  particularly  as  the  I/O  counts  on  the  die 

approach  optimum  values  for  MCM  use  (as  opposed  to  I/O  starved  PCB-optimized 
ICs),  it  will  become  increasingly  difficult  to  meet  the  requirements  for  MCM 
interconnects  using  normal  metals.  Although  appropriate  signal  propagation  has 
^  been  shown  to  be  practical  at  values  as  high  as  25-50Q,  HTSC  interconnects 

offering  values  far  below  the  desired  =  10^  range  (even  with  ultra-high 
interconnect  densities  at  linewidths  of  W  =  Tfim  or  so)  should  prove  ideal  for  high 
performance,  complex  digital  MCMs.  Further,  these  HTSC  intercoimects  should 
not  show  the  signal  anomalies  due  to  complex,  frequency  dependent  seem  at  the 
few  percent  level  in  small  copper  lines  (eg.  Figures  50, 1.17,  etc.),  which  would  make 
them  attractive  for  mixed-mode  (A/D)  applications  as  well. 
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2.1 .  Transit  Response  Characteristics  of 
Smaii  MCM  Interconnect  Lines 

2.1.1.  Background 

The  report  is  a  continuation  of  work  presented  in  the  previous  two  quarterly  reports 
and  will  continue  numbering  equations,  references,  and  figures  from  where  they  left 
ofi.  For  reference.  Equations  1-24  were  in  the  March  1992  Quarterly  Report,  along 
with  References  1-4  and  Figures  1-15.  Equations  25-45  are  in  the  June  1992 
Quarterly  Report,  along  with  References  5-8  and  Figures  16-35.  Hence  the 
equations  in  this  report  begin  with  Equation  47  and  the  figures  begin  with  Figure  36. 

Described  in  the  previous  report  was  the  development  of  a  highly  accurate  Fourier 
transform  method  of  determining  the  transient  response  of  small  lossy  lines.  The 
key  to  this  analysis  is  the  use  of  both  the  complex  propagation  constant  and  the 
complex  fi'equency-dependent  characteristic  impedance  of  the  small  diameter  lines 
(eg.  center  conductor  diameters  from  0.25  and  to  1  mil  in  small  copper  coax  lines 
have  been  analyzed).  The  principal  focus  of  this  report  will  be  on  presenting 
extensive  data  on  many  different  simulation  runs  of  this  program.  However,  as  a 
reference  point  for  the  analyses,  it  is  important  to  review  briefly  what  lossless  (ideal) 
transmission  lines  behave  like  under  similar  circumstances. 
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2.1^.  Lossless  Transmission  Line  Tiransient  Response 

Consider  the  classical  source  terminated,  open-load  transmission  line  configuration 
shown  in  Figure  36.  When  a  voltage  step  of  magnitude  is  applied  at  r  =  0 
through  a  source  resistance,  to  the  input  end  of  a  lossless  =  0)  line,  for  a 
period  of  time  of  2  (before  current  reflections  from  the  load  end  of  the  line 
arrive  back  at  the  source  end),  the  line  input  voltage,  is  given  by 

n  =  g+z  fO'-0<»<V  Eq.47 

■‘in 

The  voltage  reflection  coefficient  at  the  end  of  the  line  is  given  by  (eg.  Reference  7) 


P  = 


22  +  ^0 


Eq.48 


For  this  open-circuit  case,  =  »  and  />=+!,  which  means  that  the  load 
voltage, Kj  will  be  twice  the  voltages  on  the  line,  so  that  the  initial  output  voltage 
step  is 


2K.  Z 

^2  V 

in 


Eq.49 


Defining  the  initial  voltage  overshoot  for  this  lossless  line  case  as 


V  —V- 

Overshoot  *  — — —  for  ^  <  t  <  3^^ 


Eq.50 


in 


we  can  plot  the  expected  overshoot  versus  source  resistance  for  this  lossless  line 
case.  This  is  shown  for  Z^  =  50Q  in  Figure  37  over  a  wide  range  of  generator 
resistances  (log  scale  for  IQ  ^  ^  10,000Q),  and  in  more  detail  for  near  Z^  in 

Figure  38.  These  curves  are  very  handy,  for  example,  for  investigating  the  use  of 
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Lossless  Transmission  Line 
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self-terminating  lines  (these  are  lossy  lines  where  the  of  the  line  is  compensated 
for  by  reducing  below  Z^).  The  problem  with  this  is  that  if  the  longest  lines  are 
properly  compensated  by  reducing  R^  from  50Q  to  25Q,  for  example,  then  shorter 
lines  (which  have  small  series  resistance  and  hence  are  nearly  lossless)  can  give 
excessive  overshoot  (eg.  33.3%  for  R^  =  25Q  in  Figure  38).  Not  only  can  this 
overshoot  itself  cause  problems  with  IC  inputs,  but  large  overshoots  are 
accompanied  by  a  substantial  undershoot,  given  by 

R  Z 

“  \i  t0T3T^<,<ST^  Eq.51 

(«»  + 

This  can  seriously  compromise  noise  margin  if  the  overshoot  is  serious.  Of  course, 

the  same  is  true  of  undershoot  in  the  output  pulse  waveform,  which  occurs  for  R.^  > 

Z  in  lossless  lines. 
o 

Figures  39  through  49  illustrate  these  overshoot  and  undershoot  characteristics  for 
lossless  lines,  as  a  helpful  reference  point  for  both  design  purposes  and  for 
interpreting  the  detailed  lossy  line  computer  analysis  results  in  Section  23.  The 
generator  resistances,  R.  (called  R  in  these  particular  figures),  have  been  selected 
in  this  sequence  of  figures  to  give  overshoot  values  of  +50%,  +40%,  +30%, 
+  20%,  +10%,  +5%,  +0%,  and -5%  (5%  undershoot), -10%, -20%,  and -50%. 
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2.1.3.  Small  Lossy  Coax  Tlranslent  Response  Simulation  Examples 

As  discussed  In  the  previous  quarterly  report,  the  MCM-size  copper  interconnect 
(eg.  0.25  mil  to  1  mil  conductors)  propagation  characteristics  were  evaluated 
through  the  example  of  small  copper  coaxial  lines  with  a  lossless  dielectric  of  =  3.90 
between  the  conductor  and  shield.  As  see  in  Figure  36,  this  makes  the  velocity  of 

o 

propagation  of  a  lossless  line  equal  to  =  1.518  x  10  meters/s  or  the  unit  delay 
Ihf^  =  1673  ps/inch.  For  a  12  inch  long  interconnect,  the  nominal  delay  will  be 
^  =  2.0078  ns,  or  for  a  6"  line,  =  1.004ns,  etc.  The  nominal  characteristics 
impedance  of  the  coax  lines  is  kept  at  =  50Q  nominal  for  all  line  sizes.  This 
means  that  if  the  conductivity  of  copper  were  infinite,  that  the  line  characteristic 
impedance  would  be  =  50.00000Q  at  any  frequency.  It  is  a  reasonable 
approximation  to  the  high  frequency  Z^  of  the  real  coax  lines  where  the  skin  effect 
depth  is  small  compared  to  the  conductor  sizes.  At  lower  frequencies,  the  current 
density  spreads  deeper  into  the  conductors,  increasing  the  internal  inductance 
component  of  the  line  inductance  (Equations  5, 10 , 12, 37, 41  and  46),  which  makes 
the  a  complex  frequency  dependent  quamity.  The  fact  that  the  magnitude  of 
is  higher  at  low  frequencies  than  the  nominal  50^  value  is  why  the  R.^  »  50Q  curve 
in  Figure  50  actually  shows  an  overshoot  of  2.18%  near  ^  =  3  instead  of  the 
-0.63%  undershoot  that  would  be  expected  at  r  =  3  r  ^  from  Equat^8. 

Figures  50  through  56  show  a  composite  of  a  large  number  of  CPU  hours  of  lossy  line 
simulations  using  the  code  described  in  the  last  quarterly  report  These  simulations 
are  for  room  temperature  copper  coax,  with  center  conductor  diameters  of  1  mil,  03  mils 
and  035  mils,  lengths  of  12"  (or  6"  in  Figures  55  and  56),  and  for  a  range  of  source 
resistance  (/?^  )  values  in  the  range  of  interest  Figures  50  and  51  (expanded  time 
scale)  represent  a  typical  case  for  a  6"  MCM,  a  12"  long  line  with  a  1  mil  center 
conductor  diameter  and  having  =  12.19614Q  dc  resistance  (which  would  be 
considered  satisfactory  for  modest  logic  risetimes  in  the  HTSC/normal  metal 
interconnect  MCM  tradeoffs  of  References  5  and  6).  While  the  initial  imdershoot  at 
<  3  r  ^  extrapolated  to  /  =  r  ^  for  R^^  =  50Q  Figure  50,  as  seen  in  Figure 
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51,  the  skin  effect  loss  slows  the  leading  edge  of  the  step  response  substantially. 
While  the  center  crossing  is  delayed  less  than  0.05ns,  the  =  90%  crossing  is 
not  reached  until  about  1.9ns  after  the  nominal  pulse  arrival  time. 

Figures  52  and  53,  for  0.5  mil  center  conductor  copper  coax,  12"  long  shows  how  bad 
the  situation  gets  with  higher  resistance  lines  (R^  =  48.7846Q  there).  AWth/?^  =  50S2. 
the  50%  crossing  is  reached  in  0.22ns,  but  the  90%  level  crossing  (the  nominal 
risetime  measurement  point)  is  not  reached  until  3.25ns  after  the  nominal  pulse 
arrival  time.  As  seen  in  Figure  52,  adopting  the  "self-terminating"  approach  of 
reducing  to  20Q,  for  example,  reduces  the  delay  to  the  90%  point  to  0.5ns,  but 
causes  a  20%  overshoot  at  r  «  3  r  ^  (and  which  for  shorter  lines,  as  seen  in  Figure  37 
could  reach  a  42.8%  overshoot  if  the  R.^  values  were  made  equal  for  all  lines). 
These  issues  will  be  discussed  in  greater  detail  in  the  final  report 

If  the  conductor  size  is  scaled  down  further  to  a  0.25  mil  (635fim)  center  conductor 
diameter,  the  pulse  characteristics  for  the  12"  line,  shown  in  Figure  54  are  very  bad, 
and  cannot  be  made  reasonable  even  if  the  source  resistance,  R^  =  50Q  is  nearly 
10ns,  5  times  the  nominal  propagation  delay  of  the  line. 

Shorter  lines,  such  as  the  6"  case  shown  in  Figures  55  and  56  for  a  0.5  mil  center 
conductor  diameter,  can  tolerate  somewhat  smaller  diameters.  Here  the  nominal 
R.^  =  50Q  delay  is  about  0.85ns,  nearly  as  long  as  the  line  delay.  This  can  be 
reduced  by  reducing  R.^  ,  but  only  at  the  expense  of  overshoot.  In  this  case,  the 
R^^  value  is  nearly  25Q,  which  would  appear  to  be  near  the  allowable  limit  if 
reasonably  tight  overshoot  specifications  (eg.  5  to  10%  or  so)  are  to  be  maintained. 


Applicability  of  Superconducting  Interconnection 
Te^ology  for  Hi^  Speed  ICs  and  Systems 


DARPA  90-06 


2-5 

October  1992 


1.25 


CM 


in 

o 


o> 

in 

eo 

in 

r*. 

in 

<o 

in 

in 

U) 

in 

eo 

in 

CM 

in 

in 

o 

a 

d 

d 

(O 

d 

in 

d 

d 

CO 

d 

CM 

d 

d  9 

d 

d 

d 

d 

d 

d 

d 

d 

o 

[Al  »U|ZA  'indino  d»is  linoJio  u®do 


Time  (Relative  to  Step  Input  at  Source)  (a) 


1.9eE-09  2.00E-09  2.02E-09  2.04E-09  2.06E-09  2.0eE-09  2.10E-09  2.12E-09  2.14E-09  2.16E-09  2.18E-09  2.20E-09  2.22E-09  2.24E-09  2.26E-09  2.28E-09  2.30E-09  2.32E-09  2.34E-09  2.36E-09  2.38E-09  2.40E-09 

Tima  (Relative  to  Step  Input  at  Source)  [a] 


in 

a> 

in 

eo 

in 

in 

(D 

in 

in 

■6 

in 

CO 

in 

CM 

in  »- 

in 

o 

O) 

d 

eo 

d 

h. 

d 

<o 

d 

in 

o 

d 

eo 

d 

CM 

d 

d 

o 

d 

d 

d 

d 

d 

d 

d 

d 

o 

d 

(A]  |U|ZA  *)nd)no  d*is  tinuio  umIq 


1.5E-09  2E-09  2.5E-09  3E-09  3.SE-09  4E-09  4.SE-09  5E-09  5.5E-09  6E-09  6.5E-09  7E-09  7.5E-09  8E-09  8.5E-09  9E-09  9.5E-09  IE-08  1.05E-  1.1E-08 


cr» 


to 

CM  in 

^  to 

»-  lO 

o> 

in 

CO 

in 

in 

<o 

in 

in 

!? 

in 

<o 

in 

CM 

in 

^  in  o 

o 

OI 

^  q 

a> 

d 

CO 

d 

K 

d 

<o 

d 

in 

d 

d 

CO 

d 

CM 

d 

IT" 

d  P 

d 

o 

d 

d 

d 

d 

d 

d 

d 

o 

[A]  )U|ZA  ‘tndmo  dais  linsJio  uado 


Tims  (Relative  to  Step  Input  at  Source)  [t] 


QUARTERLY  REPORT 


2.2.  Commercial  Applications 

This  past  quarter  has  seen  a  slight  shift  in  emphasis.  The  question  we  were  asking  is 
this:  What  will  it  take  the  system  designers  to  use  conventional  room  temperature 
MCM  and  then  from  there  why  will  he  use  cryoelectronic  MCM? 

An  interesting  answer  came  out  of  this  type  of  questioning.  First,  we  find  that  IC’s 
are  optimized  for  single  packages  not  for  MCM  applications.  If  we  can  optimize  IC 
for  MCM  we  will  gain  on  at  least  two  grounds  relating  to  both  yield  and  performance 
enhancement.  First  we  can  partition  chips  to  smaller  sizes,  mount  them  on  MCM 
and  not  lose  much  if  any  in  performance,  but  gain  significantly  in  price.  Large  chips 
in  production  at  extremely  low  yields  can  be  produced  at  significantly  lower  prices  if 
they  are  partitioned  at  smaller  sizes,  since  smaller  chips  will  be  produced  at  much 
higher  yields.  The  second  point  is  that  many  logic  chips  are  now  I/O  limited,  by 
producing  chips  with  area  distributed  I/O  significant  simplification  in  design  and 
enhancement  in  performance  will  result  Therefore,  we  think  IC’s  could  lead  to 
cheaper  and  higher  performance  ^sterns. 

It  then  becomes  obvious  that  where  we  optimize  IC’s  for  cryoelectronic 
applications,  we  need  to  optimize  them  also  for  MCM’s.  This  approach  could  really 
result  in  a  major  compelling  reason  to  use  cryoelectronics. 
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Section  1.0 


1.1  Introduction 

The  focus  of  the  effort  described  in  this  report  is  updating  and  extending  the  analysis 
of  the  number  of  signal  layers  required  for  a  square  2-D  MCM  fabricated  from  10 
K-gate,  0.6  cm  square  chips,  flip-chip  mounted  with  maximal  density.  The  analysis 
follows  that  of  Ref.  7  (which  treated  the  30  K-gate,  1  cm  chip  case)  in  treating  both 
300“K  and  77“K  cases,  and  both  standard  high  speed  =  lOQ  and  compromised, 
self-terminated  R  =  25Q  cases.  It  should  be  noted  that  the  text  of  Section  1.3  is 

OUtt 

largely  reproduced  from  Ref  1,  however,  the  analytical  results  obtained  in  this  report 
are  new  and  relate  to  the  previously  mentioned  extensions.  The  analysis  is  extended 
to  include  plots  of  the  required  thickness  of  the  interconnect  mat  for  various  signal 
frequencies,  as  well  as  the  number  of  signal  layers  required  to  complete  the 
inter-chip  wiring.  For  completeness,  a  short  background  (Ref.  7)  is  included  as 
Section  12. 


1.2  Background 

Since  nearly  the  advent  of  integrated  circuit  technology,  the  vision  has  been  to 
bypass  the  great  expense  and  performance  compromises  of  multi-level  packaging  by 
implementing  entire  systems  (or  major  portions  thereof)  on  a  single  IC  wafer. 
Unfortunately,  a  number  of  technical  realities  mitigate  against  such  a  wafer-scale 
integration  approach,  at  least  for  general  purpose  digital  ^sterns.  The  principal 
problem  is  yield,  and  the  very  high  overhead  costs  of  coping  with  this  through 
redundanQT,  but  perhaps  equally  important  are  issues  of  process  and  design 
optimization.  Manufacturers  of  a  specifrc  ultra  high  density  DRAM  chip  could  not 
be  expected  to  achieve  the  same  cost,  density  and  performance  levels  if  the  function 
had  to  be  created  as  a  part  of  a  larger  general-purpose  logic  design.  These,  and 
other,  practical  considerations  have  modified  the  focus  away  from  the  single 
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semiconductor  substrate  wafer  scale  integration  vision  to  the  “hybrid  wafer  scale” 
concept  in  which  virtually  the  same  density  goals  are  achieved  in  a  multi-chip 
module  (MCM)  virtually  “tiled”  with  VLSI  die,  as  illustrated  in  Fig.  1.  Multi-chip 
modules  (MCMs)  are  an  advanced  form  of  hybrid  integrated  circuit  developed  to 
meet  the  special  needs  of  complex  high  density  digital  systems.  Specifically,  MCMs 
have  large  numbers  of  integrated  circuit  (IC)  die,  each  having  a  large  number  of  I/Os 
(signal  input/output  connections),  which  makes  the  density  of  inter-chip 
interconnections  very  high. 

It  is  now  generally  accepted  that  performance,  and  potentially  cost,  reasons,  the 
coming  wave  in  electronics  packaging  is  the  multi-chip  module  (MCM).  In  MCM 
packaging  technology  [see  Refs.  1-4],  bare  integrated  circuit  die  are  attached,  both 
physically  (chip  bonded)  and  electrically  (pad  bonded),  to  a  substrate  which  carries 
some  type  of  multi-layer  interconnect  system  to  provide  for  the  power  connections 
to  the  die  and  the  large  number  of  signal  interconnections  required  between  the 
various  chips  on  the  MCM,  and  to  the  “outside  world”.  The  MCM,  which  may 
contain  anywhere  from  a  few  IC  chips  to  over  a  hundred  in  current  MCM  practice 
(and  which  mty  be  expected  to  increase  to  many  hundreds,  as  shown  in  Fig.  1,  and 
possibly  the  low  thousands  in  the  future),  represents  the  first-level  package  for  the 
IC  die.  By  eliminating  the  requirement  for  a  large  number  of  individual  IC 
packages,  the  MCM  allows  for  major  reduction  in  size  and  weight  in  electronic 
systems  (and  potentially  reduction  of  cost). 

One  of  the  key  derivative  benefits  of  MCM  technology,  which  is  a  consequence  of 
the  reduction  in  size,  is  a  major  reduction  of  signal  propagation  delay  or  signal 
latency.  In  transmission  lines,  the  propagation  velocity  of  signal  pulses  is  limited  by 
the  speed  of  light  divided  by  the  square  root  of  the  average  dielectric  constant  of  the 
dielectric(s)  between  the  signal  line  and  ground  plane(s)  [Eq.  1].  The  total  pulse 
delay  is,  of  course,  not  less  than  the  line  length  divided  by  this  velocity.  The  very 
small  distances  between  1C  chips  in  MCM's  (as  compared  to  circuit  board 
packaging)  greatly  reduce  this  propagation  delay.  Since  in  most  digital  systems  the 
maximum  clock  frequent^  is  limited  to  the  reciprocal  of  the  sum  of  the  logic  circuit 
delay  plus  this  signal  interconnect  propagation  delay,  reduced  propagation  delays 
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translate  into  higher  system  speeds,  particularly  if  the  logic  speeds  can  be  increased 
as  well.  Of  course,  the  major  direction  of  thrust  in  the  semiconductor/IC  industry 
has  been  to  increase  chip  speeds.  This  is  rather  dramatically  illustrated  by  the  fact 
that  the  fastest  supercomputer  available  in  1980,  the  CRAY  1,  had  an  80MHz  clock 
rate  (and  $10  million  price).  In  comparison,  the  Macintosh  FX  on  which  most  of  the 
work  on  this  project  was  done  was  purchased  in  1990  for  under  one  thousandth  of 
price  of  a  CRAYl,  but  in  fact  runs  at  a  40  MHz  clock  rate,  half  that  of  a  CRAYl. 
Hence  we  see  that  today,  even  in  less  expensive  CMOS  applications,  the  reduced 
interconnect  delays  afforded  by  MCMs  are  precisely  what  digital  system  designers 
need  to  translate  the  escalating  IC  speeds  into  greatly  improved  system 
performance. 

The  simplest  concept  for  a  large  MCM  would  be  to  use  normal  semiconductor 
processing  (two  to  four  layers  of  metal  interconnects  with  typically  “l/<m  signal 
linewidths)  on  a  large  (4"  to  6"  or  larger)  wafer  to  serve  as  the  inter-chip  connection 
medium  in  the  MCM.  Unfortunately,  even  a  4'  square  MCM  would  have  up  to  8" 
long  signal  lines,  which,  assuming  a  ljum  x  0.5  fim  cross  section  and  a/>  =  1.7  x  10 
ohm  cm  metal  resistivity  (for  copper;  aluminum  would  be  worse),  the  line  resistance 
would  be  =  6900  ohms.  This  is  two  or  three  orders  of  magnitude  too  large  to 
be  useable  for  interconnects  in  a  high  speed  digital  system. 

All  signal  interconnects  in  digital  systems  are,  for  better  or  for  worse,  transmission 
lines.  Signal  lines  implemented  with  normal  metal  conductors  can  usually  be 
approximated  as  an  ideal  distributed  L-C  line  with  a  conductor  resistance  in  series 
with  each  inductor  (this  ignores  the  dielectric  losses).  In  the  ideal  case  when  the 
total  signal  line  resistance  is  very  small  in  comparison  to  the  line  impedance,  R,j^^ 
<  <  Zj,,  the  waveform  integrity  is  excellent  and  the  velocity  of  propagation,  v^,  is  just 
given  by 


Vp  =  c/Vi7 


Eq.  1 
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where  e  ^  is  the  average  dielectric  constant  in  the  line.  Since  the  line  impedances 

normally  employed  are  in  the  range  of  =  50Q  (or  within  a  factor  of  2  of  this), 
clearly  the  8"  long,  W  =  1  ^m  copper  line  =  6900Q)  case  cited  above  is  far 
from  this  ideal  case.  Note  that  for  >>  Z^,  the  propagation  approaches  a 
distributed  R-C  line  characteristic  which  exhibits  so-called  “telegrapher’s  delay”  or 
signal  “diffusion”  characteristics,  with  gross  distortion,  risetime  degradation  and 
long  signal  delays  [see  Ref.  4,  Sect.  5.2,  pp.  198-211].  While  the  exact  tolerance  to 
interconnect  line  resistance  depends  on  details  of  the  application,  in  general, 
R„^<10Q  is  quite  satisfactory,  even  for  ECL-type  (Rioad  =  Z^  =  50S2 
load-terminated)  interconnection  approaches,  with  slightly  higher  ohmic  line 
resistances  tolerable  in  source-terminated  (CMOS-like)  approaches. 

The  basic  problem  with  the  use  of  normal  metal  interconnects  in  large,  complex 
MCMs  is  that  (because  of  the  substantial  metal  resistivity)  the  conductor  areas 
required  to  keep  the  end-to-end  ohmic  line  resistances  under  lOQ  make  it  extremely 
difficult  to  achieve  the  wiring  densities  needed.  In  fact,  the  only  way  to  approach  the 
wiring  densities  needed  with  normal  metal  conductors  is  to  increase  the  number  of 
signal  interconnect  layers,  often  (as  will  be  shown  in  Section  1.3)  to  impractically 
large  numbers  for  maximal-density  complex  MCMs.  Since  a  superconductor 
behaves,  to  a  large  extent,  like  a  conductor  with  essentially  zero  resistivity,  much 
finer  conductor  linewidths  (limited  mainly  by  lithography)  could  be  used  for  such 
applications.  However,  the  fabrication  of  complex  MCMs,  even  with  the  most 
convenient  and  easy  to  use  materials,  is  challenging  from  yield  and 
manufacturability  standpoints,  and  high  temperature  superconductor  (HTSC) 
materials  are  certainly  not  convenient  and  easy  to  use  from  an  MCM  fabrication 
standpoint.  Hence  the  first  question  to  address  regarding  the  application  of  HTSC 
interconnects  to  complex  MCMs  is  "Is  this  trip  really  necessary?".  This  question  is 
addressed  in  Section  1.3,  using  quantative  analysis  tools  and  techrJques  developed 
previously  in  this  effort  and  described  in  detail  in  Ref.  1. 
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1 .3  Analysis  Approach/Nominai  300°K  Cu  to 
HTSC  Comparison  Results 

The  studies  and  analysis  discussed  above,  and  in  Refs.  1  and  7,  indicate  that  there  are 
strong  advantages  for  the  use  of  high  temperature  superconductors  for  the  signal 
interccn  lections  in  large,  very  high  density  2-dimensional  multi-chip  modules  of  the 
type  illustrated  in  Fig  1.  In  particular.  Section  3.1  of  Ref.  1  describes  in  detail  the 
technique  of  analysis  of  the  number  of  signal  intercoimect  layers  versus  MCM  size 
required  to  complete  the  inter-chip  wiring  for  maximal  density  (surface  "tiled"  with 
IC  chips  with  minimal  spacing  between  them)  square  MCMs  (eg..  Fig.  1).  The 
comparison,  between  copper  (room  temperature  in  Ref.  1  and  expanded  in  Ref.7 
and  Section  1.4  here  to  include  77“K  Cu)  and  HTSC  interconnect  materials,  is  based 
on  a  number  of  assumptions.  These  assumptions  include:  1)  chip  I/O  counts 
determined  from  Rent’s  rule  as  published  by  IBM  [Ref.  2,  page  14,  Fig  1-7]  for  high 
performance  chips  (eg.  370  I/O’s  for  a  10,000  gate  chip,  taken  to  be  0.6  cm  square), 
2)  average  wire  lengths  from  a  Donath  statistical  wire  length  model  [Ref.  3,  page  64 
or  Ref.  4,  page  430]  with  Rent’s  rule  exponent  of  p  =  2/3,  3)  a  conductor  width  of 
W,  a  thickness  of  t^= W/2  [for  fabricability]  and  a  spacing  of  3W  [for  low  crosstalk; 
see  Ref.  3,  Ch.  3,  pp.  78-124],  where  W  is  the  same  for  all  conductors  and  selected  to 
meet  line  resistance  limitations,  4)  a  40%  utilization  of  theoretical  channel  space 
by  actual  interconnects  (40%  wiring  efficiency  [eg..  Ref.  3,  page  63]),  5)  the  longest 
line  does  not  exceed  two  times  the  MCM  size  (ie.,  comer  to  corner  line  with 
"Manhattan"  wiring  with  no  "wrong  way"  routing,  6)  both  dc  and  ac  (skin-effect) 
signal  line  loss  mechanisms  are  considered  [in  Ref.  1,  just  dc  and  300MHz  ac  cases 
were  treated,  here  extended  to  a  range  of  signal  risetimes/  frequencies],  7)  a  room 
temperature  bulk  resistivity  (1.70;<Q-cm)  value  is  assumed  for  copper,  with  a  77®K 
value  4.4x  smaller,  selected  to  give  a  calculated  (from  Eq.  7)  R^^  at  f  =  IGHz  equal 
to  that  measured  in  Ref.  5,  8)  point  to  point  wiring,  making  the  number  of  wires 
equal  to  half  of  the  total  number  of  I/O  pads,  and  9)  a  maximum  allowable 
resistance  for  the  worst  case  line  of  R^^  (where  R^^  =  lOQ  is  taken  as  adequately 
small  in  comparison  to  the  characteristic  impedance  for  load-terminated  or 
source-terminated  interconnects,  and  R^  =  250  is  considered  for  either  partially 
self-terminated  or  as  a  modest  performance  compromise  in  normal 
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source-terminated  interconnects).  For  these  specific  parameters,  the  required 
number  of  signal  layers,  S,  is  given  by  (fi-om  Ref.  1,  Section  3.1,  Fig.  1-4) 

S  =  Eq.2 

where  is  the  average  wire  length  (from  Donath  model)  measured  in  units  of  the 

chip  mounting  pitch,  Xj.yp,  I  is  the  number  of  VO  pads  per  chip,  W  is  the  signal 
conductor  width,  and  the  constant  5  comes  the  product  of  0.5  (wires/pad)  times  4 
(wire  pitch/W  ratio)  divided  by  40%  (wiring  efficiency).  When  looking  at  the  plots  of 
S  vs.  MCM  size  which  follow,  it  is  desirable  to  be  able  to  identify  the  copper 
conductor  linewidth,  W,  for  any  given  point  The  corresponding  tables  of  data 
include  the  W  values  in  cm  for  the  f  =  1.0  GHz  (Tr = 400ps  risetime)  case,  but  not  for 
dc  or  other  signal  frequencies.  However,  R^  is  given  in  the  tables  and  it  does  not 
depend  on  signal  fi'equency,  etc.,  so  from  Eq.  2  the  reader  can  simply  use 

W  =  SX,,,p/5R,„I  Eq.3 

to  get  the  signal  line  width  for  any  other  case.  For  a  maximal-density  MCM 
fabricated  with  a  given  chip,  the  mounting  pitch,  X^^^,  and  of  course  the  I/O  pad 
count,  I,  will  be  constant,  and  R^^  does  not  vary  very  strongly  with  MCM  size,  so  W 
largely  tracks  the  number  of  signal  layers,  S. 

In  Ref.  1,  Section  3.1,  a  comparison  for  dc  and  300MHz  ac  (the  third  harmonic  of  a 
lOOMHz  clock,  for  example,  corresponding  to  pulse  risetimes  of  the  order  of  13ns) 
signals  of  300°K  copper  and  HTSC  interconnects  for  a  maximal-density  MCM  is  was 
described.  The  analysis  indicates  that  under  the  above  assumptions,  taking 
R„.,  =  lOQ,  a  6-inch  square  MCM  tiled  with  lOK  gate,  0.6  cm  square  chips  flip-chip 
mounted  on  a  0.635  cm  pitch  would  hold  576  chips  (24  rows  and  columns),  but  would 
require  62  layers  of  room  temperature  copper  signal  lines  to  complete  the  over 
100,000  interconnects  between  chips  (actually,  that  analysis  treated  the  ac  "skin 
effect"  line  losses  favorably  in  ignoring  the  current  crowding  and  ground  plane  loss 
contributions;  a  better  estimate  of  the  number  of  signal  layers  to  keep  the  300  MHz 
line  resistances  as  under  lOQ  is  shown  in  the  analysis  presented  later  to  be  S  =  92 
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layers  for  a  6"  maximal  density  MCM).  With  the  usual  triplate  configuration,  (one 
layer  of  X  and  one  layer  of  Y  direction  signal  ines  sandwiched  between  pairs  of 
ground  [reference]  planes)  illustrated  in  Fig.  2  for  the  case  of  S  =  6  signal 
interconnect  layers,  the  total  number  of  metal  (including  ground)  layers  would  be 
93!  While  IBM  uses  63  layers  for  their  Enterprise  System  9000  modules,  these  have 
available  a  full  set  of  "E-C"  (engineering  change)  pads  on  the  surface  with  which 
patch  and  repair  strategies  can  be  executed.  (These,  of  course,  greatly  reduce  the 
number  of  chips  that  can  be  mounted  on  the  MCM  -  eg.,  121  for  the  5"  IBM  boards). 
Creating  defect-free  boards  with  93  metal  layers  would  appear,  from  a 
manufacturing  standpoint,  highly  un-feasible,  to  say  the  least. 

With  high  temperature  superconducting  interconnects,  on  the  other  hand,  this  same 
maximal  density  6"  MCM  can  be  wired  using  only  2  levels  of  signal  interconnects 
(one  "X"  and  one  "Y").  This  S  =  2  case  is  essentially  identical  to  standard 
commercial  MCM  configurations  (eg.  the  nCHIP  process)  as  illustrated  in  Fig.  3, 
except  that  for  HTSC  all  dimensions  are  reduced  by  about  lOX  to  20X.  This  is 
possible  by  reducing  the  width  of  the  signal  lines  to  levels  (eg.,  W  of  »  1.5 
limited  by  photolithography  rather  than  line  resistance  (W  of  »  S0/<m 
for  copper).  Obviously,  even  with  only  two  signal  interconnect  layers  it  will  not  be 
easy  to  make  the  1.9  miles  of  signal  lines  on  a  6"  MCM  completely  free  of  defects, 
but  with  only  one  X  and  one  Y  layer  of  lines,  it  should  be  practical  to  develop  laser 
or  other  repair  strategies  to  handle  modest  numbers  of  defects  after  completion  of 
the  boards. 

The  attraction  of  the  HTSC  interconnect  MCM,  then,  is  that  large  "system  on  a 
substrate"  sized  maximal  density  MCMs  could  be  simply  manufactured  with  only  two 
signal  layers  (using  board  repair  strategies  to  deal  with  material  and  processing 
defects).  The  payoff  for  putting  everything  in  the  same  Ist-level  (MCM)  package  is 
avoiding  the  large  inter-package  delays  involved  in  passing  through  the  second  level 
(eg.,  circuit  board  level)  of  packaging  when  the  IC  chips  must  be  divided  among 
many  MCMs.  Since  usually  the  interconnect  delays  are  not  more  than  25%  to  40% 
of  the  cycle  time,  simply  reducing  interconnect  delays  without  corresponding 
reductions  in  the  logic  delays  in  the  ICs  would  not  buy  that  much  in  the  way  of  system 
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perfonnance.  (In  principle,  however,  sharply  reducing  the  required  total  MCM  area 
and  the  number  of  MCMs  [typically  to  only  one],  and  largely  eliminating  a  whole 
level  of  packaging,  should  be  capable  of  sharply  reducing  system  costs  [offset  to 
some  degree  by  the  cost  of  cryogenic  cooling]).  In  fact,  however,  suitably  designed 
IC  chips  have  much  shorter  logic  delays  when  operated  at  liquid  nitrogen  (77°K) 
temperature  than  when  operated  at  room  temperature  (2X  to  3X  faster  for  CMOS 
is  frequently  cited).  Hence,  by  shortening  both  the  logic  delay  (which  taken  by  itself 
would  not  be  too  effective  if  the  interconnect  delay  were  not  changed)  and  at  the 
same  time  reducing  the  interconnect  delay  (by  mounting  everything  on  a  single 
maximal  density  MCM),  adoption  of  the  HTSC  MCM  approach  could  lead  to  up  to 
3X  speed  improvements  at  the  system  level  Remembering  that  if  the  speed  of  a 
$1M  machine  is  increased  3X  it  is  worth  at  least  $3M  (sometimes  much  more  than 
that),  it  would  appear  that  there  is  plenty  of  value  margin  to  pay  for  the  cost  and 
inconvenience  of  haying  to  do  cryogenic  cooling. 


1.4  HTSC/Cu  MCM  Comparisons  vs.  Signal 

Risetime/Frequency  and  Extension  to  7TK  Cu  and 
Self-Termination  Cases  Favorable  for  Cu 

It  is  important  to  remember  that  the  reason  for  the  attractiveness  of  HTSC  MCM 
packaging  in  terms  of  system  performance,  low  cost,  small  size,  etc.,  has  nothing  to 
do  with  the  intercoimects  being  superconducting  per  se.  It  is  just  the  fact  that  the 
chips  operate  faster  at  80° K  and  that  the  use  of  HTSC  interconnects  allows  all  of 
them  to  be  placed  essentially  shoulder  to  shoulder  on  a  single  very  large,  very  high 
density  MCM  which,  because  only  2  signal  layers  are  required,  should  be  reasonably 
easy  to  manufacture,  test  and  repair.  If  by  some  means  it  were  possible  to  invent 
some  other  new  technology  or  combination  of  technologies  which  could  duplicate 
this  same  MCM  size  and  density  (and  also  operate  at  cryogenic  temperatures)  and 
also  be  manufacturable,  testable  and  repairable,  then  it  could  serve  as  well.  Of 
course,  these  alternatives  may  themselves  be  more  difficult  to  develop  than  HTSC 
MCMs,  or  have  other  disadvantages,  but  they  are  worth  considering. 
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The  first  alternative  would  be  copper  interconnects  cooled  to  80°K.  For  pulse  rise 
times  in  the  40ps  to  0.4ns  range,  the  maximum  signal  frequent^  components  are  in 
the  10’  to  10^’  Hz  range.  While  room  temperature  copper  (p  =  1.70  fi  Q-cm)  has  a 
surface  resistance,  of  8.3  mS2  at  10’  Hz  or  26  mQ  at  10“  Hz,  when  copper  wire 

is  cooled  to  77®K,  falls  to  about  3.9  mQ  at  10’  Hz  or  14  mQ  at  10“  Hz,  about  a 

factor  of  2  improvement  (these  are  measured  R^^  values  fi'om  Ref.  5,  not  simply 
scaled  firom  the  «  5  x  reduction  in  resistivity  from  300°K  to  77°K).  Assuming  that 
the  conductor  thickness,  t^,  is  well  over  twice  the  skin  depth,  to  keep  the  series 
resistance  of  a  wire  of  width  W  and  length  L  to  some  maximum  value,  R^^,  we  must 
have 


2W  Eq.  4 

This  expression  ignores  the  sidewall  contribution  to  conduction,  by  taking  as  the 
wire  periphery  2W  instead  of  the  actual  periphery  of  2(W+t^)=  3W  for  a  wire 
thickness  of  t^ = W/2.  This  is  offset  by  the  fact  that  we  have  ignored  the  effect  of  the 
nonuniform  current  distribution  on  the  wire  and  the  resistance  in  the  ground  plane 
which  increases  the  ac  resistance.  Eq.  4  is  essentially  equivalent  here  to  assuming  a 
current  form  factor  of  2/3  in  calculating  skin  effect  line  loss,  and  should  be 
reasonably  accurate.  From  Eq.  4,  the  minimum  allowable  linewidth,  necessary 
to  keep  the  resistance  of  a  line  of  length  lowe*  than  R^^j^  is  given  by 


W. 


^max  ^*urf 

~2R 


Eq.5 


For  example,  for  a  6"  MCM  with  =  12"  and  taking  R^  ~  10  Q ,  for  77*’K 
copper  would  be  =  59.4  pm  (234  mils)  at  10’  Hz  or  =  213  ^m  (8.4  mils) 
at  10“  Hz. 


To  quantitatively  evaluate  the  effect  of  cryogenic  cooling  of  copper  interconnects 
and  other  options  on  the  comparison  between  HTSC  and  Cu  MCMs,  we  will  treat 
one  specific  case  under  the  various  conditions.  Let  us  consider  a  maximal-density 
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MCM  fabricated  from  lOK-gate  VLSI  chips,  0.6  cm  square,  flip-chip  mounted  on 
0.635  cm  centers  in  a  square  array  (as  in  Fig.l),  with  370  VO  pads  per  die  (from  the 
high  performance  chip  Rent’s  rule  curve  of  Ref.2,  Pg.  14,  Fig.  1-7).  Table  1  presents 
the  calculation  results  of  the  number,  S,  of  300°K  copper  interconnect  layers 
required  to  complete  the  inter-chip  wiring  with  an  R^„=10Q  maximum  line 
resistance,  as  compared  to  the  number  of  W = 1.5m  HTSC  layers  required  as  the 
number  of  rows  and  columns  is  increased  from  NR =2  to  40  (NC=4  to  1600  total 
die),  corresponding  to  MCM  sizes  (edge  dimensions)  from  XM=0.5  inches  to  10 
inches.  Note  that  the  calculational  details  on  Tables  1-4  such  as  copper  sldn  depth, 
^iuiT  required  signal  linewidth,  W  (called  XW  in  Ikbles  1  -  4)  are  for  a  signal 
frequency  of  IGHz  (or  pulse  risetime  of  about  400ps)  only,  but  the  number  of  signal 
layers,  S,  required  for  dc  or  other  frequency/risetime  values  have  been  added  as  well. 
Use  Eq.  3  to  calculate  linewidth,  W,  values  corresponding  to  other  frequency  points 
from  the  values  of  S  at  the  desired  frequency  and  R^  (which  is  a  function  of  size 
only)  from  the  tables. 

The  calculated  number  of  30(J“K  copper  signal  layers,  S,  for  R^=  lOQ  is  plotted 
against  MCM  size,  XM  (in  inches)  for  dc  and  signal  frequency  components  of 
300MHz,  IGHz,  3GHz,  lOGHz  and  30GHz  in  Fig.  4.  Using  the  usual  signal 
risetime,  Tr,  vs  signal  bandwidth,  f,  relationship  of  f=0.4/Tr  (0.36  from  the  Fourier 
transform  of  a  Gaussian  pulse,  for  example),  these  frequencies  correspond 
respectively  to  risetimes  of  about  1.3ns  (eg.,  300®K  CMOS),  400ps  (eg.,  77®K 
CMOS),  130ps  (30(rK  GaAs  or  fast  ECL),  40ps  (eg.,  7TK  GaAs  HEMT  or  Si-Ge 
ECL)  and  13ps  (eg.,  advanced  7TK  GaAs  HBT).  Inasmuch  as  15  GHz  room 
temperature  performance  with  26  ps  to  34  ps  (unloaded/loaded)  gate  delays  in 
500-gate  gate  arrays  has  already  been  demonstrated  (using  GaAs  HBlk  with  over  55 
GHz  f,  and  f^  values),  even  the  30GHz  frequency  is  within  the  range  which  will 
have  to  be  considered  for  some  late  90’s  MCM  applications.  We  note  in  Fig.  4  that 
room  temperature  copper  is  essentially  useless  in  this  frequent  range,  requiring, 
for  =  over  22  signal  layers  just  to  interconnect  4  die  in  a  0.5"  MCM. 
Probably  the  lowest  risetime  of  broad  interest  for  mid  to  late  90’s  ICs  will  be  about 
0.4  ns,  corresponding  to  IGHz  frequencies.  At  IGHz  it  would  take  94  layers  of 
300°K  copper  signal  interconnects  to  wire  a  4"  MCM  with  256  lOK-gate  chips  on  it 
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Horizontal  is  Edge  Size  of  MCM  in  inches,  Vertical  Is  Number  of  Signal  Layers  Required  to  Complete  Inter-Chip  Wiring 

Figure  4 
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(or  142  total  metal  layers  including  grounds),  as  compared  to  2  layers  of  W  =  1.8^m 
HTSC  interconnects.  Note  that  in  a  typical  "triplate"  (as  shown  for  S  =  6  in  Fig.  2) 
configuration,  the  total  thickness  of  the  interconnect  mat  for  S  layers  of  linewidth  W 
will  be  of  the  order  of 


t„.,  =  2.25SW  Eq.6 

where  the  constant  2.25  is  appropriate  for  =  3.9  with  =  50Q  as  seen  in  Fig.  2. 

Hence,  noting  in  Table  1  for  this  S=94  case,  W=83//m,  so  the  total  wiring  mat 
thickness  would  be  an  incredible  0.7  inches  (a  thickness  of  17%  of  the  4"  size  of  the 
MCM)!  As  shown  in  the  plot  of  t^^,  versus  for  various  frequencies  in  Fig.  5,  the 
situation  is  so  far  out  of  hand  for  larger  MCM  sizes  or  higher  signal  fi-equencies  that 
the  use  of  300°K  copper  is  not  worth  discussing  there. 

As  noted  above,  ciyogenically  cooling  the  copper  can  reduce  the  resistivity  by  nearly 
5x  which  substantially  reduces  the  number  of  interconnect  layers  required  for  a 
given  MCM  size  or  signal  frequency  or  risetime.  Table  2  presents  results  of  the 
calculation  for  R^=  lOQ  for  the  maximal  density  lOK-gate  MCM  example  cited 
above,  when  the  copper  interconnects  are  cooled  to  7TK.  Note  again  that  the  77’K 
skin  depth,  and  W  values  in  Table  2  are  for  f  =  IGHz  only;  use  Eq.  7  for  R^^at 
other  frequencies,  and  Eq.  3  with  R^^  and  S  from  Ikble  2  to  get  signal  linewidth,  W, 
values  at  other  frequencies.  These  R^^  =  lOQ,  77“K  Cu  results  are  plotted  in  Fig.  6, 
and  the  wiring  mat  thicknesses  in  Fig.7,  for  comparison  to  the  HTSC  MCM.  Note 
that  cooling  the  copper  reduces  the  number  of  signal  layers  for  a  4"  MCM  supporting 
0.4ns  risetimes  (f  =  IGHz  curve  in  Fig.  6)  fi'om  S = 94  at  300®K  to  S = 44  at  77“K,  and 
reduces  the  wiring  mat  thickness  from  0.7  inches  to  a,  still  thick,  but  much  more 
reasonable  t^,  =  0.156  inches.  For  rather  slow  (probably  too  slow  for  cooled  CMOS) 

Tr  =  1.3ns  signal  risetimes  (300MHz  curve  in  Fig  6),  the  R _ » lOQ  can  be  met  even 

for  an  XM = 6"  maximal  density  MCM  with  S = 44  layers  of  77“K  Cu  (with  W = yi5fim 
lines  for  t^^,  =  2.25SW  =  0.126  inches).  While  this  77'’K  copper  represents  a 
significant  improvement  over  300°K  operation,  it  is  a  far  cry  fi’om  HTSC,  where  just 
2  signal  layers  of  W  =  1.5;im  HTSC  interconnects  would  be  adequate  for  a  6" 
maximal-density  MCM,  even  out  to  extremely  high  signal  fi'equencies.  Of  course. 
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Horizontal  is  Edge  Size  of  MCM  in  inches,  Vertical  is  Number  of  Signal  Layers  Required  to  Complete  Inter-Chip  Wiring 

Figure  6 


Tmat  for  Rmax-IOft  lOK-Qate  (370  I/O)  Flip-Chip  Max-Density  MCM  Comparison  of  HTSC  to  Copper  vs.  Frequency 


Horizontal  Is  Edge  Size  of  MCM  In  inches,  Vertical  Is  Thickness  of  Wiring  Mat  in  inches  for  Zo-50Q.  Er=3.9  Dielectric 
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while  at  least  the  wiring  mat  thicknesses  for  TT’K  Cu  are  reasonable,  fabricating  an 
S=44  signal  layer  (67  total  metal  planes  including  power/ground  planes)  with  any 
reasonable  yield  without  some  extremely  effective  repair/rework  scheme  (very 
difficult  with  67  layers!)  seems  extremely  difficult.  (If  that  is  not  obvious,  just  note 
how  complex  the  S  =  6  case  in  Fig.  2  looks,  and  imagine  what  an  order  of  magnitude 
more  layers  would  be  like!)  It  should  also  be  noted  that  no  one,  in  my  awareness,  has 
demonstrated  that  a  "conventional"  copper  polyimide,  (for  example),  MCM  with  so 
many  layers  will  be  reliable  in  temperature  cycling  to  77°K. 

The  signal  layer  requirements  and  their  dependencies  shown  above  can  be 
understood  more  clearly  with  a  bit  of  analysis.  For  a  maximum  line  length  of 
and  assuming  a  line  thickness,  t^  of  half  of  the  width,  W,  with  a  normal  metal 
resistivity, /9,  the  minimum  acceptable  width  to  achieve  an  end  to  end  dc  resistance 
of  will  be  given  by 


W 


=  V2pL 

max 


Eq.7 


We  note,  for  example,  that  the  minimum  acceptable  line  width,  W,  varies  as  the 
square  root  of  resistivity,/?,  which  is  the  same/)  dependence  as  obtained  for  the  high 
frequency  skin  effect  case  of  Eq.  5,  due  to  the  dependency  of  R^^in  Eqs.  1  and  2  on 
/),  given  by 


surf 


Eq.8 


where  n^  =  An  x  10*’  in  centimeter  units  (this  equation  ignores  the  complication  of 

anomalous  skin  effect  which  alters  the  low  temperature  somewhat).  In  either 
case,  cooling  copper  from  300°K  to  77*1C  allows  a  bit  over  2x  reduction  in  W  and 
hence  2x  reduction  in  the  number  of  signal  layers,  as  was  shown  above. 


While  the  dependencies  of  line  width  and  required  number  of  signal  layers  on  metal 
resistivity  happens  to  be  the  same  for  the  "dc"  and  "ac"  (fully  developed  skin  effect 
conduction)  cases,  this  is  not  the  case  for  their  dependencies  on  R^  or  L„,.,  which 
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tend  to  be  square  root  in  nature  for  the  "dc”  case  or  proportional  for  the  "ac"  case. 
Note  that  in  both  the  300“K  and  7T*K  plots,  in  Figs.  4  and  6,  of  the  required  number 
of  signal  layers  versus  MCM  size  for  R„^=10Q,  above  an  MCM  size  of  1.2",  the 
300MHz  ac  (skin  effect  limited)  curve  rises  above  the  "dc-only"  curve.  The  point  at 
which  these  depart  (crossover  between  dc  and  fully  developed  skin  effect  curves) 
occurs  at  a  length,  obtained  by  equating  Eqs.  2  and  3  for  W,  and  substituting  Eq. 

S^orR^^as 

^  f  Eq.  9 

where  as  usual,  f  is  the  frequency  of  the  highest  frequency  components.  here  is 

interpreted  as  the  approximate  maximum  length  for  which  a  R^  resistance  line  can 
be  considered  to  be  dc  resistance  dominated  subject  to  the  above  geometric 
assumptions.  From  Eq.  9,  for  R„.,=  lOQ  and  f  =  300MHz  (^jj=1.2566E-8 

Webers/cm  ampere),  L  should  be  6.75  cm  (2.66"),  and  since  the  longest  line  for  an 
MCM  of  size  X,^  is  L = 2Xjj,  this  corresponds  to  an  MCM  size  of  X^,  =  =  1.33" 

as  seen  in  Figs.  4  and  6.  Ikking  a  more  realistic  1=10^  Hz,  as  appropriate  for  Tr  = 
0.4ns  logic,  the  "ac"  curve  pulls  away  from  the  "dc-only"  curve  for  R„.,  =  lOQ  above 
wire  lengths  of  only  2  cm,  and  the  number  of  required  layers  for  a  6"  MCM  at  300"K 
is  S  =  178  signal  layers  (268  total)  for  R^  =  10£2 ,  or  S  =  80  signal  layers  (121  total) 
if  the  copper  temperature  is  reduced  to  77®K.  For  really  fast  logic  (eg.  GaAs  or 
Si-Ge  ECL  in  the  late  90’s),  T,^  =  40  ps  or  f  =  10*°  Hz  needs  to  be  considered. 
Here,  cooled  copper  (Table  2)  would  require  252  signal  layers  (379  total)  to  achieve 
R  =  10£2  at  lOGHz  in  the  12"  lines  for  maximal  density  6"  MCMs.  These  are,  of 
course,  hopelessly  large  numbers  requiring  unpractically  thick  "wiring  mats",  even  if 
it  were  conceivable  to  yield  them  somehow  (see  Fig.  7).  Hence,  for  higher 
performance  logic  with  higher  signal  frequencies,  copper  has  very  serious  problems 
in  attempting  to  implement  maximal  density  large  MCMs. 
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A  potential  source  of  a  significant  reduction  in  the  number  of  normal  metal  (eg., 
TT’K  copper)  signal  interconnect  layers  required  would  be  to  alter  the  chip  to  chip 
electrical  interconnection  scheme  in  such  a  way  as  to  increase  the  allowable 
maximum  line  resistance,  R  K  can  be  increased  by  some  factor,  then  for 
skin-effect  dominated  fi'equencies,  the  number  of  layers  can  be  proportionately 
reduced  (but  never  below  the  dc  resistance  limit,  of  course).  When  the  dc  resistance 
is  the  limitation,  the  number  of  signal  layers  required  will  drop  as  the  reciprocal  of 
the  square  root  of  R^.  The  value  of  =  10  Q  was  obtained,  as  described  in 
more  detail  in  Ref.  1,  Section  3.1,  as  the  point  for  =  50  Q ,  where  dc  noise  margin 
in  an  internal  -  reference,  ECL-type  (load  terminated,  as  also  used  in  GaAs) 
interconnection  approach  begins  to  be  significantly  compromised,  and  which  gives  a 
comparatively  small  loss  in  the  leading  edge  amplitudes  in  source  -  terminated,  open 
circuit  load  (CMOS  -  type)  interconnect  approaches.  For  the  latter  case,  it  is  well 
known  that  if  the  source  resistance  is  correspondingly  reduced,  a  higher  line 
resistance  can  be  tolerated,  even  up  to  R„^  =  for  R^^  =  0  if  load  capacitance 
effects  are  ignored.  For  this  ideal  case  of  R^^  =  50Q  instead  of  R^  «  10  Q  4t 
would  take  substantially  fewer  numbers  of  signal  layers  to  wire  the  MCM.  While  in 
fact  there  are  serious  problems  with  this  approach,  it  is  worth  examining  the 
implications  for  MCM  wireability  anyway.  Because  the  1C  driver  output  impedance 
tends  to  be  nonlinear,  and,  of  course,  is  non-zero,  a  realistic  expectation  might  be 
that  the  line  ohmic  resistance  might  be  able  to  be  increased,  for  Z^  -  50Q ,  to  25Q 
(with  a  combination  of  IC  driver  impedance  plus  source  termination  resistor  also 
equal  to  25S2  )  and  still  maintain  fairly  good  pulse  reflection  characteristics.  This 
R  =  2SQ  ohmic  line  resistance  should  also  be  useable  with  load  terminated  lines 

ITUUC 

if  a  differential  signal  configuration  were  used  (although  this  would  increase  by  2x 
the  number  of  wires  needed). 

This  R....= 2S£2  case  is  treated  in  Tables  3  and  4  and  Figures  8  and  10  for  S,  and  9 
and  1 1  for  t^^,.  Ikble  3  and  Figure  8  show  the  results  of  the  calculation  of  the  number 
of  room  temperature  copper  layers,  while  Table  4  and  Figure  10  show  the  results  for 
77°K  cryogenically  cooled  copper  interconnects.  From  Eq.  9,  we  calculate  for 
R^ = 25Q  a  300MHz  ac  value  for  of  16.9  cm  (6.65",  corresponding  to  an  = 
3324"  MCM  size),  or  at  a  more  practical  IGHz  signal  firequency  (0.4ns  pulse 
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Horizontai  is  Edge  Size  of  MCM  in  inches,  Vertical  is  Number  of  Signal  Layers  Required  to  Complete  Inter-Chip  Wiring 

Figure  8 
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Horizontal  is  Edge  Size  of  MCM  in  inches,  Vertical  is  Thickness  of  Wiring  Mat  in  inches  for  Zo-50Q,  Er-3.9  Dielectric 

Figure  9 
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Rmax>-25n  lOK-Gate  (370  I/O)  Flip-Chip  Max-Density  MCM  Comparison  of  HTSC  to  77"K  Copper  interconnects  vs.  Frequency 


Horizontal  is  Edge  Size  of  MCM  in  inches,  Vertical  is  Number  of  Signal  Layers  Required  to  Complete  Inter-Chip  Wiring 

Figure  10 


Tmat  for  Rmax-25fl  lOK-Gate  (370  I/O)  Flip-Chip  Max-Denslly  MCM  Comparison  of  HTSC  to  77“K  Copper  vs.  Frequency 


Horizontal  is  Edge  Size  of  MCM  in  inches.  Vertical  is  Thickness  of  Wiring  Mat  in  inches  for  Zo=50ft.  Er=3.9  Dieiectric 

Figure  11 
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risetime),  L^=5.07  cm  (or  Xj^=  1  inch  MCM  size).  Note  in  both  Figures  8  and  10 
that  the  300MHz  ac  curves  of  S  vs.  depart  from  the  dc  curves  above  Xj^,=3.3", 
and  the  IGHz  ac  curves  above  about  1"  as  calculated  from  Eq.  9.  In  terms  of  signal 
layer  requirements,  the  = 25Q,  300“K  Cu  case  plotted  in  Fig.  8  is  similar  to  the 
^max”  ^  plotted  in  Fig.  6,  and  need  not  be  discussed  further  here. 

The  most  favorable  case  for  copper  is,  of  course,  the  77®K  cooled  copper  case  with 
the  relaxed  Rn„j=25Q  (partially  self-terminated  lines)  ohmic  line  resistance 
specification.  As  seen  in  Table  4  and  Figure  10,  an  =  6  inch  maximal-density 
MCM  carrying  576  lOK-gate  flip-chip-mounted  37C-I/0  VLSI  die  having  Tr= 0.4ns 
signal  risetimes  (reasonable  for  cooled  CMOS)  would  require  over  32  layers  of  77°K 
copper  intercoimects,  or  49  total  metal  layers  (including  power/ground  planes)  in  a 
wiring  mat  about  t^,,  =  1.7tnm  (0.067")  thick  to  complete  the  1.94  miles  of  inter-chip 
wiring  on  the  MCM.  By  the  mid  to  late  90’s,  cooled  short-chaimel  CMOS  risetimes 
will  likely  be  several  times  faster  than  this,  corresponding  to  signal  frequencies  in 
the  3GHz  range.  As  seen  in  Figs.  10  and  11,  the  6"  MCM  at  this  speed  would  require 
55  cooled  copper  signal  layers  of  W=4Vim  width  to  achieve  the  R„„=250,  or 
about  84  total  metal  layers  in  a  t^^, = 5.1mm  (0.20")  thick  wiring  mat.  Hence,  even  in 
this  most  favorable  case  for  copper,  the  number  of  interconnect  layers  required 
would  appear  to  be  impractically  large  from  a  manufacturabilityyyield  standpoint. 
Further,  using  such  an  R„„=25Q  self-terminated  line  approach  would,  for  good 
signal  integrity  characteristics,  probably  require  a  number  of  different  driver  types, 
each  having  different  source  resistance  characteristics,  thereby  complicating  the 
design  task. 

As  is  pointed  out  in  Ref.l,  Section  3.1,  plain  R^  =  Z^,  open  circuit  load/source 
terminated  interconnections  can  tolerate  somewhat  over  R  .  =  lOG  line 
resistances  without  major  degradation  of  the  substantial  CMOS  noise  margins. 
From  that  standpoint,  ohmic  line  resistances  of  represent,  with  true 

open-circuit  loads,  only  a  very  modest  degradation  of  pulse  propagation 
characteristics.  However,  when  the  rather  significant  input  capacitances  which  load 
the  lines  is  considered,  pulse  edges  with  source  termination  are  degraded  relative  to 
load  termination,  and  that  is  made  even  worse  when  higher  resistance, 
self-terminated  lines  are  used.  Consider,  for  example,  a  lOpF  load  capacitance  on 


Applicability  of  Superconductbg  Interconnection 
Technology  for  Hi^  Speed  ICs  and  Systems 


DARPA  90-06 


1-15 

January  1992 


QUARTERLY  REPORT 


an  =  50  Q  tenninated  line,  which  will  have  a  550  ps  rise  time,  but  with  a  = 
50Q  source-terminated  line,  the  rise  time  is  1.1  ns.  Hence,  while  source  termination 
is  convenient  and  low  power,  it  has  high  speed  performance  limitations  relative  to 
load  termination,  which  are  only  further  aggravated  as  self-termination  is 
approached. 

A  final  potential  method  of  reducing  the  number  of  copper  layers  required  to  wire 
large  maximal  density  MCMs  is  to  relieve  the  constraint  that  all  lines  have  the  same 
width.  In  IC  design,  it  is  common  to  use  narrower  (higher  Q  /cm)  lines  for  the 
(numerous)  shorter  lines,  and  the  wider  line  widths  only  for  the  longer  lines. 
Unfortunately,  in  a  MCM,  simply  increasing  the  line  width  alone  doesn’t  help  at  all, 
since  it  correspondingly  reduces  the  line  impedance,  which  makes  the  line 
proportionately  more  sensitive  to  R„,,.  However,  if  the  dielectric  layer  thickness  is 
increased  in  the  same  proportion,  the  line  impedance  will  be  preserved  and  the 
desired  relative  line  ohmic  resistance  improved.  Hence,  it  would  be  possible  in  an 
MCM  to  create  two  or  more  classes  of  signal  interconnects,  with  wide  copper  traces 
and  thick  dielectrics  for  long  lines  and  (higher  density)  narrower  traces  with  thinner 
dielectrics  for  short  lines.  While  this  would  require  some  substantial  innovations  in 
order  for  the  CAD  tools  (particularly  MCM  routers)  to  support  such  an  interconnect 
hierarchy  (and  routing  software  is  already  a  serious  problem  for  very  complex 
boards),  it  is  certainly  theoretically  possible  to  implement  such  an  approach.  While 
its  practicability  might  be  in  question,  the  fact  that  there  are  so  many  short  wires  on 
a  MCM  makes  it  of  interest.  As  opposed  to  the  maximum  signal  line  length  of 
=  12"  on  a  6”  MCM,  as  seen  in  Tables  1  -  4,  the  average  wire  length,  R^^,  for  the 
same  case  is  only  (R^^  =  4.608  chip  pitches  of  =  0.635  cm  for  =  24)  C  =  2.93 

cm  (1.15").  Hence,  from  Eq.  5  for  10’  Hz  and  R^  =  lOD,  a  7T*K  copper  line  width 
of  W=:5.74;<m  would  be  adequate  for  half  of  the  interconnect  wires.  Since  from 
Table  2,  S = 80  layers  of  W = 59/4m  cooled  copper  interconnects  would  be  needed  for 
all  of  the  wiring,  from  Eq.  2,  using  W=5.74^m  for  short  lines  would  mean  that  half 
of  the  MCM  wires  could  be  disposed  of  in  4  layers  of  signal  interconnects,  which 
would  make  a  substantial  reduction  in  the  total  number  of  signal  layers  required. 
Unfortunately,  it  is  not  quite  that  simple.  The  longest  (L^=  12")  class  of  77“K 
copper  wires  must  be  over  10  times  wider  (59.4)Um/5.74^m)  and  the  dielectric  layers 
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10  times  thicker  than  for  short  lines.  In  general,  vias  and  other  layout  features  of  a 
size  adequate  for  thick  dielectric  layers  and  large  trace  sizes  are  not  very  compatible 
with  interconnects  an  order  of  magnitude  smaller  in  feature  sizes.  (Eg.,  via  contact 
diameters  are  often  2  or  more  times  the  linewidth,  so  for  the  long  lines  they  might  be 
12Qum  to  ISQum  or  more  in  diameter,  which  is  enormous  in  comparison  to  the 
W=5.74^m  width  of  the  short  lines  discussed  above.)  Hence,  real  design  rule 
limitations  will  play  a  major  role  in  how  effective  such  a  hierarchial  dielectric 
thickness/variable  line  width  approach  would  be  in  reducing  the  number  of  TT’K 
copper  signal  interconnect  layers  required  to  wire  large  maximal  density  MCMs.  A 
wild  guess  might  be  a  2x  reduction  over  the  constant  width  case,  but  it  certainly 
depends  on  the  details  of  the  layout  design  rules  for  the  MCM  technology. 

In  summary,  several  new  normal  metal  intercoimect  MCM  technology  extensions 
have  been  considered  which,  if  developed  and  adopted,  could  potentially  reduce 
signiHcantly  the  number  of  copper  interconnect  layers  required  to  wire  large, 
maximal  density  MCMs.  If  all  of  these  (cooling  to  77*K,  using  high  resistance 
[self-terminated]  intercoimect  lines,  and  routing  in  a  hierarchy  of  dielectric  layer 
thicknesses  and  line  widths  according  to  signal  line  length)  were  used  at  the  same 
time,  and  if  the  rise  times  of  the  logic  signal  pulses  were  fairly  slow,  then  it  might 
even  be  possible  to  reduce  the  number  of  copper  interconnect  planes  to  something 
which  might  be  manufacturable  (given  the  existence  of  very  effective  in-process 
inspection  and  repair  capabilities,  and,  of  course,  assuming  that  the  very  large  stress 
levels  that  build  up  in  organic  dielectric  MCM  systems  at  cryogenic  temperatures 
don’t  give  insurmountable  reliability  problems).  However,  if  this  can  all  be  done  in 
a  simple,  two  signal  layer  HTSC  MCM,  which  will  work  with  today’s  CAD  tools 
(routers,  etc.),  which  will  work  with  any  kind  of  signal  interconnect  termination 
scheme  (eg.,  source  or  load  terminated),  and  which  will  work  with  even  the  fastest  of 
today’s  integrated  circuit  technologies,  and  indeed  with  the  fastest  expected  in  the 
foreseeable  future,  and  do  all  of  this  without  the  performance  and  design 
compromises  discussed  above,  and  further  yet  be  potentially  cheaper,  in  the  long 
term  why  would  anyone  want  to  consider  anything  short  of  the  HTSC  MCM 
technology?  The  issue  is  not  whether  HTSC  MCMs  are  worth  doing,  it  is  whether 
we  can  in  fact  master  the  technology  to  build  them. 
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Section  2.0 

Commercialization  of  High  Temperature 
Superconducting  Technology 


2.1  Background 

Discussions  with  Tony  Vacca,  Frank  Patton,  and  colleagues  in  the  industry  and 
government  labs  have  convinced  me  that  the  program  in  HTSC  MCM  should  really 
be  viewed  from  a  much  broader  perspective.  HTSC  MCM  is  in  fact  part  of  the  larger 
field  of  cryoelectronics.  Richard  Eden  had  clearly  demonstrated  that  the 
combination  of  HTSC  MCM  and  cooled  CMOS  could  eventually  yield  an  order  of 
magnitude  increase  in  performance  over  a  corresponding  room  temperature  CMOS 
system.  V^th  the  cost  of  building  a  foundry  dramatically  increasing  (where  now  it 
would  exceed  $1B)  and  with  the  cost  of  cryogenic  decreasing,  maybe  the  time  has 
come  to  start  looking  at  cryoelectronics  as  a  real  business.  As  has  been  stressed 
before,  the  survivability  of  HTSC  companies  is  at  risk  unless  a  sizeable  market  is 
found  that  these  small  companies  can  address.  Since  these  companies  need  to  cool 
the  electronics  anyway,  it  may  make  sense  that  while  they  are  developing  the  HTSC 
technology,  they  can  start  now  looking  into  the  possibility  of  building  products  which 
are  non-superconductive.  Non-superconducting  cooled  CMOS  based  system  may 
be  a  viable  business.  These  small  companies  having  the  cryoelectronics  capabilities 
(cryogenics,  testing,  packaging,  cycling,  etc.)  could  provide  new  performances  using 
current  CMOS  technology  but  optimized  to  operate  at  liquid  nitrogen  could  make 
available  to  these  companies  a  significant  market  niche  and  render  these  HTSC 
companies  as  part  of  the  semiconductor  community.  In  January  we  held  a  work  shop 
in  San  Jose  at  N-Chip  on  this  subject  Hiere  seems  to  be  a  consensus  of  this 
thinking.  We  have  spoken  to  E-Systems  and  Cray  about  using  the  E-L  machine  as  a 
test  bed  and  to  N-Chip  about  testing  the  Ross  RISC  processor  at  LN  temperature. 
Discussions  with  semiconductor  houses  and  other  companies  producing  CMOS 
systems  will  continue  to  evaluate  the  possibility  of  converting  room  temperature 
systems  into  cooled  systems.  Commercial  and  military  applications  will  be 
considered. 
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Section  1.0 

MCM  Line  Resistance  Requirements  as 
Determined  from  Lossy  Line  Transient 
Response  Characteristics 


1.1.  Background 

As  discussed  in  several  of  the  quarterly  reports  and  in  the  final  report  to  the  first 
portion  of  this  effort,  the  primary  interest  in  HTSC  interconnects  in  MCMs  stems 
from  the  difficulty  of  obtaining  adequately  low  signal  line  ohmic  resistances  using 
fineline  metal  interconnects.  In  the  initial  analysis  of  the  required  number  of  signal 
layers  versus  MCM  size  for  a  maximal  density,  Rent’s  rule  I/O  count  limited  MCM, 
a  maximum  dc  end-to-end  line  resistance  of  10  Q  was  assumed  to  give  adequate 
digital  pulse  propagation  characteristics  (dynamic  noise  margins)  for  either 
source-terminated,  open-load  (CMOS-type)  or  for  load-terminated  (ECL  -  type) 
signal  interconnect  schemes  for  assuming  transmission  line  characteristic 
impedances  of  the  order  of  =  50  Q.  It  was  subsequently  noted  that  because  of 
the  substantially  faster  clock  rates  and  signal  edge  rates  possible  with  advanced 
integrated  circuit  technologies,  such  as  short  channel  CMOS,  particularly  when 
operated  at  cryogenic  temperatures,  the  relevant  high  frequency  ac  skin  depth,  <5  can 
be  considerably  smaller  than  the  signal  conductor  dimensions,  making  the  ac 
resistance  much  larger  than  the  dc  resistance.  In  the  final  report  and  the  last 
quarterly  report  (Q4  1991),  the  plots  of  the  required  number  of  signal  layers  for 
copper  interconnects  were  made  for  dc  and  for  ac  frequencies  of  300MHz,  IGHz, 
3GHz,  lOGHz  and  30GHz,  for  total  worst  case  dc  and  ac  resistances  less  than  either 
lOG  or  25G  for  either  300“K  copper  or  77®K  copper  interconnects.  A  "rule  of 
thumb"  f  =  0.4/t,  relationship  between  digital  signal  risetime  and  equivalent  ac 
signal  frequency  was  noted  for  each  of  the  frequency  curves  in  the  four  figures. 
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While  this  simple  association  of  risetime,  signal  frequency  and  a  maximum  allowable 
ac  resistance  specification  for  acceptable  signal  quality  is  a  convenient  extension  of 
the  simple  maximiim  allowable  dc  trace  resistance  concept,  it  is  not  as  rigorous  or 
quantitative  as  the  dc  specification.  The  problem  is  that  while  a  pulse  train  at  some 
clock  frequency  f^  with  a  risetime,  t^  certainly  has  strong  frequency  components  up 
to  at  least  3fp  or  Sf^,  or  up  to  frequencies  of  the  order  of  f  =  0.4/tp  these  higher 
frequency  components  do  not  appear  with  100%  of  the  signal  amplitude,  and  hence 
a  given  level  of  attenuation  of  these  components  is  not  necessarily  as  serious  as  an 
attenuation  of  the  fundamental,  or  dc  component,  by  the  same  amount.  Hence,  a 
specification  that  the  ac  skin  effect  line  resistance  at  the  third  or  f^fth  harmonic  of  f^, 
or  at  f  =  0.4/t^  must  be,  for  example,  ^  lOQ  may  be  unecessarily  severe. 

The  only  way  to  quantatively  resolve  the  issue  of  what  conductor  size  is  required  for 
a  given  length  of  interconnect  in  order  to  achieve  adequate  dynamic  noise  margin 
characteristics  under  specified  signal  risetime,  etc.,  conditions  is  to  know  accurately 
the  time-domain  (as  opposed  to  frequency  domain)  response  of  the  lossy  line.  For 
the  case  of  linear  driver  and  load  characteristics  (ie.,  driver  and  load  impedances 
which  are  linear  in  signal  amplitude  and  not  time  dependent,  but  which  may  be 
frequenty  dependent,  and  impulse  response  for  the  lossy  line  path  may  be  derived, 
and  the  output  waveform  obtained  as  the  convolution  of  the  impulse  response  with 
the  input  waveform.  In  the  spirit  of  the  rest  of  the  analysis  work  on  the  HTSC  MCM 
advantages/requirements,  we  would  like  to  come  up  with  simple  analytic  expressions 
or  standard  numerical  calculation  methods  that  anyone  can  verify,  extend  or  apply  to 
other  cases  of  interest  "in  the  comfort  of  their  own  home  or  office”,  so  to  speak.  This 
is  as  opposed  to  simply  buying  the  appropriate  CAE  tools  for  a  suitable  engineering 
workstation  and  running  these  to  get  the  transient  lossy  line  outputs  for  specific 
cases  of  interest.  (Quad  Design,  of  course,  has  only  encouragement  for  those  with 
the  budget  resources  and  inclination  to  take  the  easy  way  out  and  use  their  XTK 
[Crosstalk  Tool  Kit]  suite  of  CAE  software  tools  including  the  XFX  [Crosstalk  Field 
Extractor]  and  XNS  [Crosstalk  Network  Simulation]  tools  to  take  this  approach.) 
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This  bit  of  "tongue  in  cheek"  aside,  it  should  be  noted  that  the  purpose  of  this  effort 
is  to  derive  insight  which  guides  the  design  process  by  fitting  a  limited  number  of 
accurately  numerically  treated  specific  cases  and/or  analytic  expressions.  These  case 
may  be  idealized  and  simplified  as  required  to  illuminate  the  essential  elements  of 
the  design  problem.  This  is  quite  different  fi-om  the  (far  more  difficult)  problem  for 
which  the  Quad  Design  tools  were  developed.  These  tools  were  designed  to  execute 
very  computationally  intensive  design  evaluation  tasks  within  an  acceptably  short 
time  on  reasonably-priced  workstations.  For  example,  they  can  handle  multiple 
conductors  of  arbitrary  geometry,  multiple  segments  of  transmission  line  of  different 
geometry  in  complex  paths,  and  (the  very  difficult  case  of)  real  nonlinear  IC  output 
driver  and  input  receiver  lossy  transmission  line  sources  and  loads.  They  are 
designed  to  execute  these  analysis  in  an  "automatic"  fashion  from  board  or  MCM 
design  databases  on  large  numbers  of  real,  complex  interconnect  nets  and  furnish 
the  analysis  results  (eg.  time  domain  response  data,  delay  times,  rise/fall  times, 
undershoot/ringing  characteristics,  etc.)  for  use  in  higher  levels  of  design  evaluation 
(eg.,  timing  verification  or  logic  simulation).  This  requires  an  extremely  efficient 
implementation  of  computational  approaches  in  order  to  give  engineering  quality 
results  in  an  acceptable  time  without  requiring  a  CRAY  C-90  supercomputer  to  run 
the  job.  This  necessitates,  of  course,  both  some  cleverness  in  obtaining  very  high 
calculational  efficiency,  and  some  engineering  approximations  in  favor  of  faster 
simulation  time  where  the  sacrifices  in  aoruracy  are  not  generally  significant  for  real 
designs.  However,  for  an  analytical  fitting  of  the  type  we  are  interested  in  here, 
focus  may  be  extended  into  such  "extreme"  regions  in  which  the  engineering  tools 
may  not  furnish  the  desired  levels  of  precision.  For  this  reason,  a  much  slower,  but 
"proof  quality"  result  for  simple  single  lines  with  linear  source  and  load  is  better  for 
our  purposes. 
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1.2.  Simple  Approaches  to  Lossy  Line  Transient  Response 
1.2.1.  Coaxial  Line  Example 

Figure  1  illustrates  the  methods  used  for  calculating  the  resistive  line  loss 
contribution  to  transmission  line  response  for  the  simple  case  of  a  coaxial  line.  For 
the  purposes  of  our  analysis,  we  will  assume  the  dielectric  within  the  line  is  ideal, 
that  is,  lossless  and  non-dispersive  (a  frequency-independent  dielectric  constant,  Sj), 

with  a  value  of  e^.  =  3.90  taken  for  the  numerical  example  (eg.,  Si02  or  typical 

polyimide).  If  the  coaxial  center  conductor  (signal  line)  and  outer  shield  (ground) 
were  ideal,  that  is,  had  zero  resistivity,/?  =  0,  then  the  line  would  have  a  distributed 
external  inductance  (inductance  per  meter)  of 


(Rb 


ford  =  0 


valid  where  the  skin  depth,  d,  given  by  (Ref.  1) 


d  = 


\<ofiJ2p 


is  very  small  compared  to  the  dimensions  of  the  conductor  (R,  and  R^  -  R^,  in  Fig. 
1).  It  is  assumed  here  that  the  dielectric  and  conductor  materials  are  non-magnetic, 
that  is/i  =  =  4jr  X  10”^  henrys/meter.  For  this  case,  the  capacitance  per  meter 

of  the  coaxial  line  is  given  by 


2;re  e 
In/R/Rj 
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where  =  8.85418782  x  10"^^  farads/meter,  and  the  high  frequency  (<5»»0) 

characteristic  impedance  of  the  coaxial  line  with  center  conductor  radius  of  R,  and 
inner  shield  radius  of  Rj,  is  given  by 

Z„  ^  ^/DC  ^  ^In  y  (4) 

'  V 

where  rj  =  =  376.7303134  ohms  is  the  impedance  of  free  space.  For  the 

case  treated  numerically  here,  with  Rj,  =  6.59200577  x  10“^  meters  gives,  with  Ej. 
=  3.90,  =  50Q. 

1.2.2.  Rdc  -  Only  (Low  Frequency)  Lossy  Coax  Response 

At  the  other  extreme  to  the  d  »  0  case  discussed  above  is  the  low  frequency  "dc"  case 
where  the  skin  effect  penetration  depth,  d,  is  large  in  comparison  to  the  conductor 
dimensions  (  d>Rg  and  d  >Rq  -  Rj,).  Here  the  signal  current  density  is  not 

confined  just  to  the  thin  surface  regions  of  the  conductors  (darkly  shaded  in  Fig.  T), 
but  rather  the  ourent  flows  uniformly  over  the  conductor  area.  In  this  case  there  is 
additional  line  inductance,  due  to  the  presence  of  magnetic  fields  within  the 
conductors.  For  the  roimd  center  conductor  in  the  coax  this  "internal"  inductance  is 
given  by 

i0Td»Ra  (5) 

Also,  the  effective  radius  of  the  outer  conductor  used  to  calculate  the  external 
inductance  in  Eq.  1  should  be  about  *•  ^R|^  +  R^j  /  2.  The  total  low 

frequency  inductance  is  just  the  sum  of  from  Eq.  1  and  from  Eq.  5.  For  the 
line  dimensions  cited  above,  with  R^  =  8.128  x  10-5  meters  outer  shield  radius 

(Fig.  1),  the  coaxial  line  inductance  would  increase  from  =  329368  nH/meter 
at  high  frequency  (<5  0)  to  Ldc  =  401.409  nH/meter  (50nH/meter  increase  from 
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the  center  conductor  and  the  remainder  due  to  the  increase  in  effective  shield 
radius,  R^,')  at  low  frequency.  This  would  increase  the  line  impedance  from 

Zq  =  50Q  at  high  frequency  to  55.198Q  at  low  frequency  (a  10.4%  increase).  The 
velocity  of  propagation,  given  for  the  lossless  case  by 

Vp  =  1  /  VLC  =  c/VF  (6) 

(where  c  =  2.99792458  x  10^  meters/second  is  the  speed  of  light)  will  also  be  altered 
by  the  change  in  inductance,  from  lArp  =  6.587363ns/meter  (167319ps/inch)  at  high 

frequency  to  about  1/Vp  =  7.272ns/m  (184.7ps/inch)  nominally  at  low  frequency  if 

Eq.  6  were  still  valid  (with  L  =  +  Lj^)  for  this  d  >  >  and  S  >  >  R^  -  Rj,  case. 

However,  in  the  presence  of  substantial  line  loss  the  pulse  propagation  gives  rise  to 
distortion  and  attenuation  of  input  pulses  which  make  it  more  difficult,  in  general,  to 
define  a  velocity  of  pulse  propagation  for  the  line.  The  purpose  of  the  lossy  line 
response  analysis  is  to  calculate  the  actual  transient  output  waveform,  which 
includes  both  velocity  and  dispersive  attenuation  distortion  effects  on  pulse 
waveforms. 

As  was  discnissed  in  the  first  quarterly  report  on  this  contract  (reproduced  as 
Appendix  1  of  the  final  report  of  10/91),  the  dc  line  resistance  is  of  obvious 
importance  for  the  case  of  load-terminated  interconnect  lines  due  to  gross  loss  of  dc 
noise  margin  because  of  the  ohmic  line  resistance/load  termination  resistance 
voltage  divider  action.  K  the  dc  line  resistance  exceeds  10  or  15  £2,  the  ECL  noise 
margins  are  seriously  compromised  (unless  full  differential  signal  transmission  is 
used).  As  was  discussed  in  that  report,  while  there  is  obviously  no  problem  with 
ohmic  line  resistance  at  dc  if  the  CMOS  like  source-terminated/open  load 
configuration  is  used,  the  tolerance  for  line  resistance  for  ac  or  fast  transient  signals 
is  only  a  little  better  than  for  the  load-terminated  case. 

Because  of  skin  effect,  the  ac  line  resistance,  as  illustrated  in  Fig.  1,  is  generally 
much  higher  than  at  dc.  At  first  glance,  this  could  lead  one  to  think  that,  for  fast 
risetime  signals  (which  carry  very  high  frequency  components),  the  dc  line 
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resistance,  Rjp,  would  not  be  important.  However,  when  it  is  noted  that  for  a  1  foot 
length  (30.48  cm)  of  coax,  the  current  pulse  to  a  step  input  on  the  tp^  =  2ns  length 
of  line  is  A  t  =  2tp^  =  4ns  wide,  we  see  that  even  if  the  input  pulse  risetime  is  40ps, 
representing  frequency  components  near  lOGHz  (where  the  skin  effect  depth  in 
room  temperature  copper  is  only  d  =  0.656i«m),  the  full  envelope  of  the  4ns  current 
pulse  requires  a  low  frequency  bandpass,  for  negligible  (eg.  2.5%)  droop,  of  IMHz 
(assuming  single  pole  rolloff),  at  which  frequency  d  =  65.^m  (2.58  mils;  larger  that 
the  conductor  dimensions  in  the  coax  example). 

The  effect  of  (only)  on  the  transient  response  of  the  Imil  diameter  center 
conductor  =  50Q  coax  was  evaluated  for  lengths  of  6",  12",  18"  and  24"  in  the 

50Q  source  resistance  open  circuit  load  configuration  using  PSpice.  For  purposes  of 
calculation  the  internal  inductance  effects  were  ignored  (Z^  kept  at  50Q  and  Vp 

unchanged).  The  transmission  line  was  broken  into  lOps  long  segments  (eg.  200 
segments  to  model  the  12"  line,  400  for  the  24",  etc.)  with  a  R^g  =  0.0546034Q 
resistance  placed  in  series  with  each  ZO  =  50,  TD  =  lOP  transmission  line  segment. 
PSpice  transient  simulations  with  a  lOps  risetime  input  step  were  carried  out  using  a 
5ps  or  less  iteration  step  size.  Run  times  were  slow  (»  100  seconds/ns  for  the  12" 
line  on  a  MacIIfx),  as  expected  for  SPICE,  in  comparison  to  the  streamlined  Quad 
Design  tools  which  also  include  skin  effect  and  are  therefore  more  accurate  as  well. 
These  PSpice  R^g  (only)  lossy  line  response  simulations  are  plotted,  along  with 
other  simulations,  in  Figs.  7  through  15.  Since  they  ignore  skin  effect,  the  R  jg  curve 
is  generally  the  highest  curve  in  each  of  the  figures.  Fig.  12,  for  a  24"  (60.96cm)  line 
length  of  the  Imil  copper  coax,  illustrates  the  Rjg  (only)  response.  The  initial 
output  step  amplitude  after  the  tpg[  =  4ns  line  delay  is  8038%  of  the  input  voltage 
step  (it  would  be  100%  if  the  line  were  lossless,  of  course).  This  initial  step 
amplitude  is  given  from  the  total  dc  ohmic  line  resistance,  R^,  (eg.  R^g  = 
21.84 136Q  for  the  24"  example  of  Figs.  12  and  13)  by 

=  v^e  for  R,,  (only)  (7) 


Applicability  of  Superconducting  Interconnection 
Te^ology  for  Hi^  Speed  IC’s  and  Systems 


DARPA  90-06 


1-7 

March  1992 


QUARTERLY  REPORT 


where  is  the  input  step  amplitude  and  tp^  is  the  propagation  delay  time  of  the 
line  (4ns  in  Fig.  12).  As  seen  in  Fig.  12,  the  step  response  for  (only)  increases 
approximately  linearly  (it  is  probably  exponential-like)  over  the  next  2  tpj  to  a 
second  break  point  at  t  =  tp^  (from  the  initial  input  step)  at  which  the  output 
amplitude  is  given,  at  least  approximately,  by 


/  V  -012 

2 

fv  (t  )) 

out  '  Dd  ' 

V  .f30  =  V.  e 

z 

o 

V. 

out\  pd/  in 

V  J 

1  “  / 

for  (only)  (8) 


where  the  quantity  (tp  j)  in  the  exponential  expression  is  found  from  Eq.  7. 

The  numerical  value  of  0.12  in  Eq.8  was  determined  by  fitting  the  data  for  the  four 
line  lengths  analyzed,  and  gives  an  accuracy  for  VQyj(3tpj)  of  ±0.005%  for  the  cases 
treated.  Beyond  this  t  =  3  tp^j  break  point  (at  12ns  in  Fig.  12),  extrapolates 
more  or  less  linearly  to  essentially  at  t  =  5  tp^  (20ns  in  Fig.  12). 

Eqs.  7  and  8  for  the  R^j^  (only)  response  break  points  apply  most  accurately  to  the 
cases  of  most  interest;  that  is,  where  R^^  us  substantially  smaller  than  (so  the  line 
transient  response  is  favorable).  Just  for  interest,  a  high  resistance  case 
corresponding  to  12"  of  about  0.3  mil  coax,  with  a  total  Rjp  =  lOO.OQ  or  Rjc/Zp  = 
2  for  a  12"  line  was  analyzed  and  found  to  fit  Eqs.  7  and  8  fairly  accurately  (025%  for 
Eq.7  and  1.9%  error  for  Eq.  8)  even  though  it  is  in  the  highly  lossy  (erfc-like  step 
response)  region  (very  slow  risetime  beyond  the  initial  step).  As  will  be  discussed 
later,  the  combined  skin  effect  plus  dc  ohmic  resistance  lossy  line  response  can  be  at 
least  approximated  by  a  convolution  of  the  skin  effect  response  (with  R^^.  =  0)  with 
this  Rjg  (only)  response.  While  the  R^j^  (only)  curves  in  Figs.  7-15  were  obtained 
using  PSpice  simulation,  simple  linear  interpolation  between  the  breakpoint  values 
given  above  (at  t  =  tp^  in  Eq.  7,  t  =  3tpj  in  Eq.  8  and  •»  Vj^  at  t  =  5tpj)  should 
be  quite  adequate  to  define  the  R (only)  response  without  PSpice. 
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1.2.3.  Skin  Effect  Only  (Rdc  =  0)  Lossy  Coax  Response 

In  general,  due  to  the  non-uniformity  of  current  density  on  the  surface  of  conductors 
in  transmission  line  configurations,  2-d  solver  numerical  methods  are  required  to  get 
accurate  values  for  the  frequency-dependent  ac  (skin-effect)  resistances  of  lines  (ie, 
to  determine  the  current  "filling  factor").  For  the  particular  case  of  a  coaxial  line, 
illustrated  in  Fig.  1,  the  magnetic  field  strengths  are  uniform  over  the  surfaces  of  the 
center  conductor  and  the  inside  of  the  outer  shield,  and  the  field  strengths  are  easily 
calculated  (falling  as  the  inverse  of  the  radial  distance).  Hence,  as  shown  in  Fig.  1, 
the  frequency-dependent  ac  (skin-effect)  resistance  can  be  easily  calculated  for  this 
case,  assuming  the  frequency  to  be  high  enough  so  that  the  penetration  depth,  d  in 
Eq.2,  is  much  smaller  than  the  conductor  dimensions  (fully  developed  skin  effect 
condition).  The  frequency-dependent  skin  effect  resistance,  in  Fig.  1,  is  given  by 

a «  R..  (R.-R.)  (9) 

(for  Rjjp  in  Q/meter),  and  the  internal  inductance,  Lj^,  due  to  the  skin  effect  is  given 
for  the  angular  frequency,  (o,  by  [Ref,  1] 

Lin.  =  R.C  for  <5  <  Ra  .  (Ro  -  Rb)  (10) 


(for  Ljjjj  in  henrys/meter,  but  the  inverse  Fourier  transform  necessary  to  convert  the 
frequency  domain  results  to  a  time  domain  step  response  is  known  in  closed  form  for 
this  fully  developed  skin  effect  case.  The  distributed  inductance,  of  a  lossless 
line  gives  rise  to  the  distributed  reactance  of  the  line 

(lossless  line)  (11) 
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For  the  lossy  case,  this  becomes 


Z  =KuL.+  kwL.  +  R 

ser  J  ext  ■'  int  a 

Z  =  j<U  L  ,  +  R 

scr  J  ext  s  X" 


(full  skin  effect)  (12) 


where  the  R^  Vjw  term  includes  both  the  real  part,  R^^  from  Eq.  9  and  the  inductive 
part,  Lj^p  from  Eq.  10.  Note  that,  from  the  45®  angle  of  this  skin  effect  term. 


Rs  = 


VTr. 


This  definition  of  R^  is  the  same  as  that  used  by  Rubin  in  Ref.  3.  Gardiol  (Ref.2) 

does  an  analysis  for  the  case  of  an  air  coax  line  (instead  of  filled  with  a  dielectric 
constant  ),  and  unfortunately  uses  the  symbol  in  place  of  fj  (see  Eq.  4). 

Converting  his  notation  to  the  notation  used  here  and  in  Ref.  3,  Gardiol’s  distance 
parameter,  z,  equals  x  Vi^,  where  x  is  the  physical  distance  down  the  line  in  meters, 

and  Gardiol’s  constant  Ic"  (as  in  the  term  -kz  Vja  in  his  Eq.  1025)  is  equal  to 
Rj  /  (2  Zq  Vc~  ),  where  is  the  line  impedance  (eg.  500).  Gardiol  gives  the  step 

response  for  this  fully  developed  skin  effect  case  at  position,  z,  and  time,  t,  as 


u  (z,t)  =  erfc 


fort  s  z/c  (14) 


where  erfc  is  the  complementary  error  function,  plotted  in  Fig.  2  for  convenient 
reference,  and  c  is  the  speed  of  light  Converting  this  step  response  expression  to  our 
notation. 


u  (x,t)  =  erfc 


4Z 

O 


fort&x/c  „  (15) 
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where  x  is  the  physical  distance  from  the  source  end  of  the  line  in  meters  and  c^^  = 
c  is  the  speed  of  light  in  the  dielectric.  is  given  for  the  coax  case  by  Eqs.  13 
and  9  as 


R.  K 


(16) 


for  <5  <  <  Rg  abd  d  <<  R^  -  Rj, ;  Eq.  15  may  also  be  written  using  a  time 
normalizing  parameter,  t^^j,  as 


u(x,t)  =  erfc 


for  t'  >  0  (17) 


where  t'  =  t  -  x/c^ff  is  the  time  measured  relative  to  the  nominal  pulse  arrival 
time,  and  t^^j  must  be  given  by  (from  Eq.  15) 


VI 


(18) 


The  complementary  error  function,  plotted  in  Fig.  2  for  reference,  may  be 
approximated  for  small  values  of  the  argument  as 


erfc  (x)  «  1  -  1.12837878X  for  x  <<  1  (19) 


Eq.  19  is  0.1%  accurate  to  x  =  0.054, 0.5%  to  x  =  0.123, 1%  to  x  =  0.174, 2%  to  x 
=  0.246, 5%  to  X  =  0.396  and  10%  to  x  =  0.574. 


Since  most  of  the  electrical  response  issues  concern  the  "tail"  of  the  response 
function  where  erfc(x)  is  nearing  unity  (x  is  small),  Eq.  19  is  a  very  useful 
approximating  function  for  cases  of  interest.  (Remember  that  in  Eq.  17  the  square 
root  of  time  is  in  the  denominator,  so  large  time  is  small  x.)  In  this  "tail"  region  of 
Fig,  2,  the  complimentary  error  function  may  be  approximated  by 
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ertc(x)- 0.48235  (20) 


This  approximation  is  good  to  ±1.25%  for  1.596  ^  x  ^  4.332,  ±5%  for  1.30 
^  X  :S  4.704,  and  ±  10%  for  1.05  ^  x  <  5.02.  Since  the  lossy  line  step  function 
expression  in  Eqs.  14, 15  and  17  involve  a  square  root,  it  would  be  more  convenient 
if  the  exponent  of  x  in  Eq. 

20  were  2,  instead  of  the  more  accurate  1.9.  A  reasonably  good  fit  for  this  tail  region 
of  the  erfc  curve  is  given  by 

erfc  (x)  »  0.419  ^  for  1  <  x  s  3.2  (21) 


This  expression  is  accurate  to  ±  10%  for  1.035  s  x  ^  3.174,  or  +10%,  -20%  for 
0.7292  ^  X  ^  3.478.  Eq.  21  is  very  useful  for  the  initial  rise  portion  of  the  fully 
developed  skin  effect  response  function .  For  example,  using  the  approximation  of 
Eq.  21,  the  voltage  response  of  Eq.  17  becomes  simply 

u(x,  t)  »  0.419  )  for  0  <  t '  t^^i  (22) 


where  t'  =  t  -  tpj  =  t  -  x  Vij.  /  c  and  tj^i  is  given  by  Eq.  18.  Eq.  22  is  very  useful  for 

picking  off  the  10%  or  20%  response  points  for  a  10%  -  90%  or  20%  -  80%  risetime 
calculation,  for  example.  The  other  end  of  the  curve,  where  the  90%  or  80%  points 
would  be  picked  off  can  be  approximated  by  using  Eq.  19  in  Eq.  17,  or 

u  (x,t)  -  1  -  1.12838  for  t '  >  33 1^^,  (23) 


which  is  1%  accurate  for  t'  >  33  t^^j  (near  the  80%  point),  and  0.25%  accurate  near 
the  90%  point  (t'>  126  t^gj).  For  convenience  in  other  regions  where  these 
approximations  are  not  very  accurate,  the  erfc  (1/Vt'  j  )  function  is  plotted  versus 
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t'/tggj  for  the  ranges  of  0  s  t'/t^gi  ^  35  in  Fig.  3,  0  ^  ^  ^ 

0st'/tsei^750inFig.  5. 

Standard  inexpensive  engineering/scientific  mathematical  packages  such  as 
Bimillenium  Corporation’s  "HiQ"  or  Wolfram  Research’s  "Mathematica"  for 
Macintosh  computers  have  efficient  built-in  erfc  (x)  functions,  as  well  as  solvers  that 
can  provide  inverse  values  (ie.  extract  crossing  times).  (The  material  for  Figs.  2-15 
here  were  calculated  using  "HiQ",  except  for  the  R  [only]  curves  [from  PSpice]  and 
Quad  XNS  simulations.) 

1.2.4.  Comparison  of  Idealized  Coax  Cases  to  XFX  and  XNS  Results 

Getting  all  of  the  stray  ^  (or  or  larger)  factors  in  Fig.  1  or  the  foregoing 

equations  (particularly  Eqs.  13  -  18)  exactly  right  is  fairly  demanding,  so  it  is  very 
valuable  to  check  these  expressions  again'  standard  engineering  tools  for 
correctness.  Since  the  coax  case  shown  in  Fig.  1  and  applicable  to  Eqs.  1, 3, 4, 5, 9 
and  16  (the  other  equations  are  general)  is  the  only  one  not  requiring  2-d  solver 
numerical  techniques  to  get  R^^  or  R^  values  accurately  (Eqs.  9  and  16),  this  case  was 
used.  The  first  step  was  to  extract  the  R^^  and  R^  (ac)  lossy  line  parameters  (along 
with  characteristic  impedance  etc.)  using  the  XFX  part  of  the  Quad  Design  XTK 
tool  kit.  Note  that  a  coax  line  is  a  veiy  difficult  case  to  solve  with  high  accuracy  on  a 
rectangular  conductor  grid  solver  (which  is  optimal  for  "real"  MCM  or  circuit  board 
interconnects  like  stripline  or  microstrip).  The  materials  values  taken  were  p  — 
1.70  X  10“®S2-m  for  300®K  copper  and  Cj.  =  3.90  for  the  dielectric,  with  the  = 
SOD  line  geometry  scaled  from  a  1  mil  diameter  center  conductor  (Ra  =  1.27  x 
10“^  m,  Rj,  =  6.592  x  10”^  m  and  R^  =  8.128  x  10”^  m,  as  defined  in  Fig.  1). 
Overall  grid  sizes  for  the  XFX  calculation  of  200  x  200, 400  x  400, 600  x  600  and  800 
X  800  points  were  evaluated  for  mesh  sensitivity.  The  R^  and  results  were  quite 
accurate  from  XFX;  eg.  Z^  =  49.98Q  even  for  the  coarsest  (200  x  2(X))  mesh  case. 
Rjg  values  for  the  4  mesh  sizes  were  0.35829Q/cm,  0.35836D/cm,  0.35847Q/cm,  and 
03585 ID/cm  respectively,  versus  a  value  of  0.35943245D/cm  calculated  analytically 
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(Fig.  1),  corresponding  to  absolute  errors  from  0.26%  for  the  800  x  800  mesh  up  to 
030%  for  the  200  x  200  mesh.  The  numerical  problem  is  a  lot  tougher  for  the  ac  skin 
effect  resistance,  reported  as  R^jj  by  the  XFX  program,  because  it  is  "texture" 
sensitive  and  the  texture  of  a  small  circle  is  poorly  represented  by  a  rectangular  grid, 
the  R.;;  values  obtained  varied  from  0.44001  for  the  200  x  200  grid,  0.44442  for  400  x 
400, 0.49029  for  600  x  600  and  0.55318  for  the  800  x  800  grid,  with  the  units  indicated 
as  QVns/cm.  These  values  were  (finally)  compared  to  a  nominal  analytic  value  for 
R^.  =  0.48847978  for  this  quantity,  so  that  the  four  numerical  results  represent 
errors  of  -10%,  -9%,  +0.37%  and  +13.25%  respectively  (not  unreasonable, 
considering  the  highly  unfavorable  geometry  for  the  numerical  approach). 


A  more  difficult  problem  was  understanding  just  what  XFX  was  reporting  for  this 
Rgij  quantity  and  how  it  relates  to  R^  or  R^^^  (Eqs.  16, 13  and  9).  For  this  example, 
we  have,  from  Eq.  16,  Rj  =  2.184548  x  10“^  Q  Vi/mctcr.  Unfortunately,  despite 
the  similarity  of  notation  to  R^  in  Ref.  3  or  in  Fig.  1  or  Eqs.  13  or  16  here,  the  XFX 
definition  of  R^jj  is  given  by 


(XFX) 


Vox  10"’ 


V2xl0 


(24) 


Hence,  the  Eq.  16  value  of  R^  =  2.184548  xlO  ^  Q  V  ns/m  divided  by  V2  x  10“’ 
becomes  R^jj  =  48.8479783  Q  V  ns/m  or  the  R^ij  =  0.4884798  Q  V  ns/cm  value 
cited  above  (which  is  in  the  center  of  the  range  of  the  XFX  calculated  values). 
Because  the  XNS  simulator  accepts  as  input  the  line  parameters  extracted  by  XFX, 
and  we  wanted  to  use  analytically  acciurate  values  for  R^  in  the  comparisons  between 
the  transient  waveform  calculations  here  and  the  XNS  waveform  simulations,  the 
analytic  R^jj  value  was  used  as  input  to  XNS  for  all  simulations  of  the  1  mil  center 
conductor  copper  coax  line  (Figs.  6-15). 

As  a  check  of  the  fully  developed  skin  effect  (only;  Rj^.  =  0)  step  response 
expressions  (eg.,  Eqs.  14-18),  XNS  was  run  with  R^jj  =  0.48848  and  R^j^  =  0 
(obviously  non-physical  unless  the  center  conductor  of  the  coax  had  a  room 
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temperature  superconducting  core)  for  a  12"  (x  =  03048m)  line  length  of  = 
50f2  coax  with  a  50f2  source  resistance  and  open  load,  using  a  Ips  risetime  input  and 
a  Ips  step  size.  From  Eq.  18,  t^gj  =  11.0839ps  for  this  line  length,  and  the  step 
response  expression  in  Eq.  17  was  evaluated  with  a  Ips  step  size  for  comparison  to 
the  XNS  simulation.  The  results  are  compared  in  Fig.  6.  Aside  from  an  apparent 
artifact  in  XNS  that  the  first  8  time  points  (past  tpj  =  2ns,  of  course)  were  printed 
as  zero  values  (they  all  had  values  under  10%  of  the  input  step  and  were  somehow 
ignored  in  the  output),  the  erfc  (l/Vt'/t^gj)  and  XNS  values  compared  very 

accurately  to  20  or  30  ps,  with  XNS  falling  below  the  analytical  expression  by  2.7% 
to  2.9%  at  larger  times.  The  reason  for  this  small  difference  has  not  been  identified, 
but  it  could  involve  either  computational  shortcuts  in  XNS  or  the  fact  that  XNS 
treats  the  full  (incident  and  backscattered)  propagating  wave  conditions,  while  Eq. 
17  assumes  only  a  forward  propagating  waveform.  In  any  event,  the  agreement  in 
Fig.  6  is  quite  good  between  XNS  and  the  analytic  functional  form  of  Eq.  17. 

1.2.5.  Lossy  Line  ITansient  Response  Calculations  with  Both  Rdc  and  Skin 
Effect  -  Practicai  Approach 

As  was  discussed  in  section  1.2.2,  because  even  fast  risetime  step  signal  require  low 
frequency  components  to  accurately  represent  them  (frequencies  for  which  the  skin 
effect  depths  will  typically  exceed  the  conductor  dimensions),  both  the  dc  (R^c)  and 
the  skin  effect  (R^)  losses  must  be  included  to  obtain  accurate  results  for  real  lossy 
lines.  Unfortunately,  the  inverse  Fourier  transform  to  convert  lossy  line  frequency 
response  to  time  domain  transient  response,  which  is  known  in  closed  analytic  form 
for  the  fully  developed  skin  effect  case  (Eqs.  14-18),  is  not  known  in  closed  analytic 
form  for  the  mixed  R^g  +  skin  effect  case.  In  fact,  in  general  even  the  frequency 
domain  response  is  not  well  known  in  the  "crossover"  region  between  low 
frequencies  (where  R^  dominates)  and  high  frequencies  (where  R^  dominates), 
(ie.,  where  the  skin  effect  depth  [d  in  Eq.  2]  is  comparable  to  the  conductor 
dimensions).  The  frequency  dependence  is  known  analytically  for  the  round  wire 
center  conductor  in  coax,  but  the  Bessel  function  solutions  [Ref.  1  or  Ref.  2]  are 
fairly  complicated.  In  Section  2.4  of  the  October  1991  Final  Report  for  this  contract. 
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approximate  analytic  expression  for  the  frequency  dependence  for  microstrip  and 
stripline  conductor  geometries  were  developed  which  would  be  very  helpful  for  use 
with  numerical  methods  for  obtaining  the  inverse  Fourier  transform  (IFF)  to  get  the 
time  domain  response.  Analytic  closed  form  expression  for  the  IFT  step  response 
are  unlikely,  however. 

In  the  real  world  of  engineering  difficult  theoretical  or  computational  problems 
must  somehow  be  dealt  with  in  a  "good  enough"  practical  way  that  gives  adequate 
accuracy  for  the  intended  engineering  purpose  of  the  tools.  The  approach  used  in 
the  Quad  XNS  tool  to  simultaneously  include  both  skin  effect  and  line  loss 
contributions  is  alluded  to  in  Ref.  3.  Figs.  7-9  compare  the  results  of  the  XNS 
calculation  for  the  12"  open  coax  example  described  previously  (done  with  a  lOps 
input  risetime  and  lOps  step  size)  with  the  fully  developed  skin  effect  expression 
from  Eq.  17  [erfc  (l/Vt'/t^gj)],  and  with  the  R^g  (only)  results  from  the  PSpice 

simulations  discussed  in  Section  122.  As  would  be  expected,  both  the  R^g  (only) 
and  skin  effect  (only)  curves  are  too  optimistic,  since  each  omits  an  important  part 
of  the  loss  (see  Fig.  9,  for  example). 

The  real  line  response  would  involve  some  type  of  convolution  between  the  R^g 
(only)  and  erfc  (l/Vt'/t^gj)  skin  effect  (only)  curves  in  Fig.  9.  Since  the  Rjg  curve  is 

fairly  flat,  a  good  first  guess  at  the  combined  response  is  just  to  do  a  point  by  point 
multiplication  of  the  R^jg  (only)  curve  with  the  erfc  (l/Vt'/t^gj)  curve  (Eq.  17).  This 

combined  curve,  noted  as  R^jg'erfc  in  Figs  7-15,  falls  in  fact  very  near  the  result  from 
the  Quad  XNS  tool.  In  Fig.  9,  for  example,  they  practically  superimpose,  at  least  out 
to  about  t  =  3.3ns  (t'  =  1.3ns).  Beyond  t  =  7.7nsec  the  XNS  result  exceeds  even 
the  skin  effect  (only)  curve.  However,  note  that  this  is  in  the  "far  tail"  region  where 
the  response  is  beyond  95%  of  its  final  value,  which  usually  is  of  little  interest  in 
digital  systems. 

Figs.  10  and  11  show  the  results  for  an  op)en  circuit  18"  piece  of  the  same  1  mil  center 
conductor  copper  coax.  Again,  the  R^jg^erfc  curves  from  this  analysis  fall  virtually 
on  top  of  the  Quad  XNS  curves  out  to  t  =  6ns  and  the  error  does  not  exceed  about 
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\25%  out  to  20ns.  There  is  also  a  small  («*  23ps)  shift  of  part  of  the  initial  rise  (eg., 
at  the  50%  crossing  point)  between  the  XNS  and  R^*erfc  curves,  but  after  tp^  = 
3ns  this  is  not  a  significant  error. 

Figs.  12  and  13  show  similar  good  agreement  for  the  24"  coax  case,  particularly  out 
to  t  =  93Sns.  Again,  however,  the  maximum  error  (aside  from  the  small,  ***  33ps, 
leading  edge  time  difference  )  is  only  about  135%.  For  the  shorter,  6"  open  coax 
case  shown  in  Figs.  14  and  15,  the  agreement  is  excellent  between  our  R^*erfc 
curve  and  XNS  out  to  t  =  1.35ns,  but  even  then  approximation  is  not  bad  out  to  t  = 
3ns  or  so.  Beyond  this,  XNS  is  substantially  too  high,  but  since  the  response  is  over 
95%  here,  the  error  is  not  important  for  digital  applications. 

1.2.6.  Comment  on  Results 

The  goal  of  being  able  to  calculate  reasonably  accurate  lossy  line  step  response 
curves  for  the  simple  coax  line  case  with  only  standard  inexpensive  computational 
tools  has  been  achieved.  The  product  of  the  fully  developed  slrin  effect  (Eq.  17) 
curve  times  the  R^j^  (only)  curve  (which  can  be  adequately  defined  without  PSpice 
simulation  for  the  open  circuit  case  by  interpolating  between  the  break  points  given 
in  Eqs.  7  and  8)  gives  reasonable  -  lookmg  results  which  agree  very  well  with  the 
XNS  simulations.  Probably  the  only  way  to  get  more  accurate  "proof  quality"  line 
transient  response  results  is  to  carry  out  a  frequency  domain  calculation  as 
accurately  as  possible  and  then  use  numerical  inverse  Fourier  transform  methods  to 
get  the  time  domain  results.  Computing  hardware  can  do  Fourier  transforms  so 
quickly  now  that  this  may  be  an  attractive  approach  for  fully  linear  cases,  however  it 
does  not  generalize  to  the  non-linear  characteristics  of  many  drivers  and  receivers  of 
practical  interest. 
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To  briefly  touch  on  the  signiflcance  of  the  calculational  results,  rather  than  the 
methods,  we  can  see  in  Figs.  7-lS  that  even  for  the  open  line  case,  line  dc  resistances 
of  the  order  of  lOQ  are  important  to  avoid  gross  compromise  in  step  response.  We 
will  be  looking  into  the  risetime  and  dynamic  noise  margin  implications  of  these 
results  in  the  future. 
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Section  2.0 

Commercialization  of  High  Temperature 
Superconducting  Technology 


2.1.  Background 

During  this  quarter,  emphasis  was  put  on  cryoelectronic  technology  market 
development  issues.  In  earlier  reports  we  discussed  the  significance  of  this 
technology  to  the  commercial  and  defense  sectors.  What  we  attempted  to  do  in  this 
quarter  is  to  better  understand  the  cryoinfrastructure.  Specifically,  what 
product/service  would  a  cryoelectronic  company  provide  to  survive.  It  is  clear  that 
the  semiconductor  houses  recognize  the  performance  increase  gained  by  CMOS  as 
it  is  chilled.  It  is  also  recognized  by  system  houses  that  chilled  CMOS  system  could 
yield  factor  of  two  or  higher  in  system  performance.  So  the  question  is  what  would 
it  take  to  get  this  technology  to  be  accepted  by  the  industry.  Recent  cost  reduction 
and  performance  improvement  in  cryogenics  indicate  that  this  part  of  the 
cryoelectronics  may  not  be  the  long  pole  in  the  tent.  Furthermore,  the  next 
generation  of  CMOS  foundries  indeed  vdll  cost  nearly  one  billion  dollars.  So  it 
would  seem  a  major  near  term  leap  frog  can  take  place  if  the  various  pieces  of  the 
technology  can  be  brought  together  and  particularly  if  the  cryopart  can  be 
transparent  to  the  system  designer.  Obviously,  it  is  difficult  to  educate,  at  least  in 
the  near  term,  mai^  system  houses  on  how  to  use  cryogenics.  However,  there  are 
companies  such  as  Sn  and  Conductus  which  alreatfy  have  ciyolectronic  know  how, 
due  to  their  experience  in  superconductivity.  These  companies  however  do  not  have 
much  experience  in  the  digital  industry. 
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Our  job  has  been  to  work  with  them  to  help  orient  to  the  right  market  niche.  It 
seems  that  there  is  room  for  a  business  which  supplies  a  system  solution  to  those 
system  houses  intending  to  use  cryoelectronics.  That  includes  device  optimization, 
CAE-tools,  packaging,  integration  and  assembly  as  well  as  testing.  Another  business 
niche  may  be  to  work  with  workstation  houses  to  take  current  and  next  generation 
products  and  optimize  for  ciyoapplication  using  some  architecture.  Although  only  a 
few  days  have  been  spent  in  this  area,  we  believe  significantly  greater  time  is  needed 
with  these  people  to  get  the  IC  industry  moving.  Potentially  new  generations  of 
workstations  can  be  realized  within  a  relatively  short  term  without  having  to  wait  for 
the  emergence  of  next  generation  CMOS. 
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Section  1.0 


Transient  Response  Characteristics  of  Smaii  MCM 
Interconnect  Lines,  Including  Complex,  Frequency 
Dependent  Characteristic  Impedance  Effects 


1.1.  Background 

As  discussed  in  the  background  for  the  last  quarterly  report,  the  critical  driving  force 
toward  the  use  of  high  temperature  superconducting  (HTSC)  interconnect  lines  in 
large,  maximal  density  multi-chip  modules  (MCMs)  is  line  resistance.  Adequately 
small  end  to  end  ohmic  line  resistances  (much  smaller  than  the  characteristic 
impedance,  Z^)  are  only  achievable  with  normal  metal  interconnects  by  making  the 
conductors  quite  large  (~lmil  for  room  temperature  copper).  These  large 
conductor  diameters  (or  widths  for  stripline  or  microstrip  conductors)  severely  limit 
the  number  of  interconnections  which  can  be  completed  on  each  layer,  and  hence 
require  excessive  and  impractical  numbers  of  signal  layers  to  implement  large, 
maximal  density  MCMs  (see  references  5  and  6). 

The  impetus  to  develop  HTSC  MCM  technology,  then,  is  a  fairly  quantitative  issue, 
depending  quite  specifically  on  just  what  conductor  dc  and  ac  (skin  effect)  resistance 
levels  are  tolerable  in  digital  MCM  interconnect  applications.  The  last  quarterly 
report  (Eqs.  1-24,  Figs.  1-15,  References  1-4  and  associated  text.  Sections  1.1-1Z.6) 
addressed  analytical  approximations  which  allow  estimates  of  the  transient  response 
of  normal  metal  interconnect  lines  to  be  generated  with  minimal  computational 
effort,  and  compared  the  results  for  the  specific  case  of  small  coaxial  lines  to  those 
generated  by  the  Quad  Design  XTK  suite  of  commercial  CAE  tools.  It  was  noted  in 
that  report  that  both  the  analytic  methods  and  the  commercial  CAE  tools  involve 
significant  approximations  (in  order  to  keep  run  times  acceptable  in  an  engineering 
environment  for  the  commercial  tools),  which  have  noticeable  impact  on  the 
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accuracy  of  the  results.  In  particular,  substantial  differences  were  seen  between  the 
methods  explored  in  the  "tail"  region  of  the  rise  characteristic,  particularly  above 
90%  to  95%  of  the  final  value,  even  though  the  basic  assumptions  behind  the 
calculations  were  the  same.  Fortunately,  the  20%  to  80%  or  10%  to  90%  regions  of 
the  signal  waveforms  are  of  the  greatest  interest  in  most  digital  MCM  or  PCB 
applications,  so  this  is  not  a  serious  concern  in  a  commercial  CAE  environment. 
However,  for  mixed  mode  MCM  applications,  more  detailed  transmission  line 
characterization  could  be  required,  and  in  any  event,  it  would  be  comforting  to  know 
what  really  is  the  true  lossy  line  transient  response,  in  order  to  evaluate  the  impact 
of  the  various  approximations  necessary  to  facilitate  the  "practical"  solutions 
discussed  in  the  previous  report. 

This  report,  then,  simply  picks  up  where  the  previous  report  stops.  Hence,  in  order 
to  avoid  endless  usage  of  "equation  or  reference  number  (x)  of  the  previous  report", 
we  will  simply  continue  the  numbering  of  equations  (starting  here  at  Eq.  25),  Figures 
(starting  here  at  Fig.  16),  and  References  (starting  here  at  Ref.  5)  from  where  the 
last  report  stopped. 
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1 .2.  Non  -  Approximate  Forms  of  Basic  Lossy  Transmission 
Line  Equations 

The  basic  approach  to  the  lossy  transmission  line  problem  discussed  in  the  last 
quarterly  report  (Eqs.  1-24,  Figs.  1-15  and  Refs.  1-4)  is  the  perturbation  method  (see 
section  4.6  of  Ref.  2),  with  the  assumption  of  small  losses  such  that,  for  lossless 
dielectrics  (only  conductor  losses),  the  complex  propagation  constant,  is  given  by 


y  +  a  +jp  «  jo)  VET  +  ^ 


Eq.  25 


and  characteristic  impedance,  Z.,  by 


Eq.  26 


These  are  approximations  to  the  general  expressions  for  lossless  dielectrics/lossy 
conductor  lines  for  the  complex  propagation  constant 


y-^JjQ)C(R+  j(oL) 


Eq.  27 


and  the  complex  characteristic  impedance 


^0  = 


R  +  icdL 


'  R+joL 
j(oC 


Eq.  28 


For  completeness  (for  the  benefit  of  those  for  whom  a  few  years  may  have  passed 
since  their  last  transmission  line  course),  the  frequency  dependent,  complex 
propagation  constant,  y ,  is  used  in  the  voltage  and  current  propagation  equations 


V(z)  =  V^*e-^^+V-e- 


Eq.  29 
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=  /+e"’'*  +  I- 


nz)-i:e 


Eq.  30 


where  the  +  and  -  superscripts  correspond  to  propagation  in  the  +x  and  -x 
directions,  and  the  complex,  frequency  dependent  characteristic  impedance,  Z^,  is 
defined  by 


K  -K 

<7  _  ®  ® 

o  o 


Eq.  31 


As  noted  in  section  1.2.3  of  the  last  report  (Equations  9-23),  the  low  loss 
assumptions  of  Eq.  25  and  26,  taken  along  with  the  assumption  of  pure  skin  effect 
loss  (R^  =  0),  allows  the  frequency  domain  signal  propagation  expression  to  be 
analytically  Fourier  transformed  back  to  the  time  domain,  giving  the  erfc  function 
step  response  solution  (Eqs.  14-19).  Analytic  Fourier  transforms  are  not  available  to 
convert  frequency  -  domain  results  calculated  using  the  general  lossy  line  equations 
(Eq.  27  and  28),  particularly  when  both  dc  and  ac  (skin  effect)  line  resistance 
contribution  are  included.  Hence,  the  strategy  utilized  here  is  to  avoid  the  use  of 
approximations,  solving  in  the  frequen<y  domain  the  general  transmission  line 
equations  with  complex,  frequen<7  dependant  y  and  and  then  using  numerical 
techniques  to  transform  the  results  back  into  the  time  domain.  While  this  technique 
can  always  be  used  to  model  the  linear  elements  (successions  of  transmission  line 
segments,  boundary  shunt  impedances,  etc)  in  a  problem,  we  will  go  ahead  and 
assume  that  the  voltage  generator  (F)^)  source  impedance,  and  the  Z,g,  load 
impedance,  Z2,  are  also  linear.  This  allows  us  to  solve  the  entire  problem  in  the 
frequency  domain  and  do  a  single  inverse  Fourier  transform  (iFFT)  to  get  the  result 
(as  opposed  to  XNS,  which  allows  nonlinear  drivers  and  loads  by  using  the  impulse 
response  of  the  line  segments  in  a  discrete  time  simulation  of  the  circuit). 
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The  problem  solved  is  illustrated  in  Fig.  16.  Since  almost  all  of  the  variables  are 
complex  frequency  dependent  values  (eg.,y ,  Z^,  Zj,  Zy„ ,  and  all  V’s  and  I’s),  for 

simplicity  in  notation  they  will  not  be  specially  marked  to  show  this.  (It  is  fairly  clear 
that  the  quantities  such  as  R’s,  Us,  C‘s  and  distances  such  as  jc  or  /  are  scalars). 
Treated  is  a  voltage  generator,  with  source  impedance,  Z,„,  connected  through  a 
single  lossy  transmission  line  segment  of  characteristic  impedance  Z^  and  length  1,  to 
a  load,  Z2.  The  unknowns  are  the  voltages  at  both  ends  of  the  line,  at  the  source 
end  and  V2  at  the  load  (frequently  we  are  most  interested  in  the  load  voltage,  V^,  but 
this  necessitates  first  finding  V^).  Numerically,  of  course,  the  solution  is  carried  out 
for  all  frequencies  before  carrying  out  the  iFFT  to  get  the  time  domain  solution. 

The  approach  to  the  frequency  domain  solution  is  illustrated  in  Fig.  16.  First,  the 
load  impedance,  Z2,  is  transformed  back  to  an  equivalent  impedance,  Z^,  located  at 
the  source  side  of  the  line,  using  the  general  transmission  line  equations  (Eqs. 
27-31).  A  bit  of  messy  complex  arithmetic  is  saved  by  the  fact  that  these 
transformations  were  done  may  decades  ago  and  are  published  in  many  standard 
handbooks  (see,  for  example,  Ref.  7).  This  Z2  to  Zj  impedance  transform  through 
length  X  is  just 


Zj  Zj  cosh  (yx)  +  Z^  sinh  (yx) 
cosh  (yx)  +  Z^  sinh  (yr) 


Eq.  32 


where  for  our  case,  x  =  1.  One  of  the  cases  of  interest  is  for  an  open  circuit  line 
(Z2  =  «  ),  in  which  case  Eq.  32  becomes 

for  Z2  open  Eq.  33 


—  =  coth()a) 


Applicability  of  Superconducting  Interconnection 
Technology  for  Hi^  Speed  IC’s  and  Systems 


DARPA  90-06 


1-5 

June  1992 


QUARTERLY  REPORT 


Having  transformed  the  load  impedance,  Z^,  back  to  an  equivalent  at  the  source, 
we  can  use  the  simple  voltage  divider  equation  to  obtain  the  voltage  at  the  source 
end  of  the  line  as 


Eq.  34 


or,  if  Zj  happens  to  be  infinite,  such  as  for  an  open  circuit  line  at  dc  (f  = 

0).  To  get  the  voltage  at  the  load  end  of  the  line,  K2,  having  calculated  the  voltage  at 
the  source  end,  V^,  from  Eq.  34,  we  simply  use  the  transform 


cosh  (yx)  +  (Z^/Zj)  sinh  (yx) 


Eq.  35 


We  can  use  the  iFFT  operation  to  go  from  the  frequency  domain  V2  to  obtain  the 
time  domain  output,  V2(t). 

It  should  be  noted  at  this  point  that  while  the  frequency  domain  solutions  based  on 
Eqs.  27,  28  and  32-35  are  deceivingly  simple  and  concise,  they  involve  complex 
functions  (sinh,  cosh,  coth,  sqrt,  etc.)  of  complex  arguments,  which  would  be  a  bit  of 
a  chore  to  deal  with  if  the  programming  language  doesn’t  support  complex  variables 
and  have  the  requisite  functions  in  its  library.  For  this  work,  the  program  "HiQ" 
from  Bimillennium  Corp.,  101  Albright  Way,  Los  Gatos,  CA  95030,  (800)  488-8662, 
running  on  either  a  Mackintosh  Ilfic  or  a  Powerbook  170  was  used  for  calculations 
and  most  graphics  output  Presumably,  Mathematica  would  also  work  for  this. 


I 
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1.3.  Coaxial  Line  Frequency  Dependence  of  Resistance  and 
Inductance  -  Kelvin  Function  Solutions 

As  was  discussed  in  section  122,  because  even  fast  rise  time  step  signal  require  low 
frequency  components  to  accurately  represent  them  (frequencies  for  which  the  skin 
effect  depths  will  typically  exceed  the  conductor  dimensions),  both  the  dc  (J(^)  and 
the  skin  effect  (R^)  losses  must  be  included  to  obtain  accurate  results  for  real  lossy 
lines.  This  creates  one  problem  in  that  the  inverse  Fourier  transform  to  convert 
lossy  line  frequency  response  to  time  domain  transient  response,  which  is  known  in 
closed  analytic  form  for  the  fiilly  developed  skin  effect  case  (Eqs.  14-18),  is  not 
known  in  closed  analytic  form  for  the  mixed  R^  +  skin  effect  case.  We  are  getting 
around  that  problem  by  using  numerical  methods  for  taking  the  iFFT  to  obtain  the 
time  domain  response.  Unfortunately,  in  general,  even  the  frequency  domain 
response  is  not  well  known  in  the  "crossover"  region  between  low  frequencies 
(where  Rj^  dominates)  and  high  frequenci<^  (where  R^  dominates),  (ie.,  where  the 
skin  effect  depth  [d  in  Eq.  2]  is  comparable  to  the  conductor  dimensions).  In  Seaion 
2.4  of  Ref.  5,  approximate  analytic  expression  for  the  frequency  dependence  for 
microstrip  and  stripline  conductor  geometries  were  developed  which  would  be  very 
helpful  for  use  with  numerical  methods  for  obtaining  the  inverse  Fourier  transform 
(IFF)  to  get  the  time  domain  response  for  these  conductor  geometries. 

Fortunately,  the  frequency  dependence  of  the  ac  resistance,  and  internal 
inductance,  L^f ,  of  the  round  wire  center  conductor  in  coax  are  known  analytically 
as  a  function  of  frequency  for  frequencies  from  dc,  where  Eq.  5  is  valid,  up  to  the 
region  where  the  high  frequency  approximation  of  Fig.  1  or  Eqs.  9  and  10  becomes 
valid  (eq.,  the  skin  depth  d  [Eq.  2]  becomes  much  less  than  the  center  conductor 
radius,  R^).  In  the  example  of  a  SOD  coax  with  =  3.90  dielectric  constant,  the 

shield  inner  radius,  R^  is  S.190SS572  times  the  center  conductor  radius,  and  hence 
the  shield  contributions  to  R^  and  are  only  16%  of  the  total.  Hence,  it  will  be 
fairly  accurate  if  we  simply  scale  the  results  for  the  center  conduaor  by  the  factor 
[1  +  RJRi,],  that  is: 
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and 


t^ac^coax) 


/{^(center  conductor 


Eq.  36 


(coax) 


^int  conductor 


Eq.  37 


As  noted  in  Fig.  1,  the  dc  resistance  (per  meter)  of  the  center  conductor  is  given  by 


^dc  )  Eq.  38 

and  its  dc  internal  inductance  by  Eq.  5.  The  frequency  dependence  of  the  internal 
inductance  and  ac  resistance  of  the  center  conductor  is  carried  out  in  terms  of  the 
dimensionless  ratio,  between  the  center  conductor  radius,  ,  and  the  (frequency 
dependant)  skin  depth,  d ,  as  (Ref.  2) 


t  =  VI  R/d  Eq.  39 

As  described  in  Section  4.53  of  Ref.  2  (Gardial),  the  real  and  imaginary  parts  of  the 
center  conductor  impedance  per  unit  length  are  given 


\v  [bei  an + v>er  a)r 


and 


^  ter  (?)1 

VW  (6ei'(0r  +  Iter'(S)] 
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where  Lj^ ^  I  Sn  as  given  in  Eq.  5.  The  functions  ber{ i;)  and  h«(C)  are  the 

Kelvin  functions  (Ref.  8),  which  are  the  real  and  imaginary  parts  of  complex  Bessel 
functions.  In  the  HiQ  program  these  zero  order  Kelvin  functions  are  called 
respectively  as  ber  (0,zeta)  and  bet  (0,zeta).  The  primes  indicate  first  derivatives 
with  respect  to  C  •  which  numerically  was  executed  with  a  3-point  derivative  with  a 
step  size  of  0.02. 

Fig.  17  shows  a  plot  of  (coax)  and  (coax)  for  the  1  mil  center  conductor, 

=  3.90,  nominal  =  50Q  coax  [with  room  temperature  copper 

(p  =  1.70;<Q  /cm)  center  conductor  and  shield],  example,  calculated  using  Eqs. 
36-41.  Note  that  in  the  "skin  effect  only"  (R^  =  0)  approximation  of  Eqs.  9  and  10, 
the  Rgg  curve  goes  to  zero  atf  =  0  and  approaches  infinity  as  the  frequency,/, 
approaches  zero.  The  Kelvin  function  calculated  curves  in  Fig.  17  appear  to  exhibit 
the  proper  behavior  and  values. 

One  practical  problem  in  using  the  Kelvin  function  solutions  is  that  these  complex 
functions  end  up  being  evaluated  eight  times  (4  ber  and  4  bei)  for  each  frequency 
point  (to  give  the  functional  values  and  3-point  derivative  values),  which  is  a  great 
deal  more  computational  effort  than  the  simple  approximate  expressions  of  Eqs.  9 
and  10.  Further,  at  very  high  frequencies  (where  (;  can  get  very  large),  the  numerical 
acctiracy  of  the  functional  evaluations  could  become  questionable  (particularly  the 
derivatives).  To  solve  both  problems,  simple  analytical  approximations  to  these  ber 
Ibei  expressions  were  developed.  The  function  to  be  approximated  for  R^ 
evaluation  purposes,  RacBer  is  given  by 


RacBer  (  0  “ 


Eq.  42 
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and  the  function  for  evaluation  is 


LintBer{^)  = 


lbei’a)f  +  lber'a)f 


Eq.  43 


Both  of  these  functions  approach  the  value  of  1/V2  for  very  large  values  of  ^  Figure 
18  shows  the  basis  for  the  simple  approximation  for  RacBer  ( for  ^  large: 


^  1  .  0.5011275  ^ 

RacBer  ( 0  ss  —  +  (g_o.475)~  ^ 

This  approximation  is  accurate  to  about  1  part  in  10^  for  ([  >  10  ,  and  0.02%  for 
(;  >  5.  The  approximation  for  LintBer  5  is  illustrated  in  Fig.  19: 


LintBer  ( C)  =  “R - 0-5156  ^  ^ 

^  (C- 0.677)^ 

This  approximation  is  accurate  to  about  2  parts  in  10^  for  C  >  10  and  about  0.02%  for 
5  >  6.  Fig.  20  compares  the  actual  RacBer  and  LintBer  functions  of  Eq.  42  and  43 
(points)  with  these  approximate  forms  of  Eqs.  44  and  45  (solid  curves). 
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1.4.  Frequency  Dependence  of  (Complex)  Characteristic 
Impedance  of  Small  Copper  Coax  Lines 

The  coaxial  line  capacitance  per  meter,  C  (Eq.  3),  and  inductance  external  to  the 
conductors  themselves,  (Eq.  1)  are  used  to  calculate  the  nominal  characteristic 
impedance,  (nominal)  of  the  line  (Eq.  4)  from  =  ylL^^/C  .  At  high 

frequencies,  where  the  skin  depth  is  very  small  in  comparison  to  the  conductor 
dimensions,  the  actual  characteristic  impedance  approaches  this  nominal  value.  In 
general,  however,  the  inductance  is  larger  thanL^  by  the  amount  of  the  inductance 
due  to  magnetic  fields  within  the  conductors  themselves,  that  is, 

L  =  L^  +  L^^  Eq.  46 

In  the  numerical  frequency  domain  calculations,  the  value  from  Eqs.  41  and  37 
(using  Eq.  45  for  $  >  10 )  is  added  to  (Eq.  1  to  give  the  total  L  value  to  go  into 
Eqs.  28  (for  Z^)  and  27  (for  y).  Similarly,  the  value  from  Eq.  40  (with  Eq.  44  for 
10)  is  used  in  Eq.  36  to  include  the  shield  loss,  and  that  R^^  (coax)  value  goes  in 
for  R  in  Eqs.  27  and  28,  along  with  the  C  value  from  Eq.  3  to  calculate  the  complex, 
frequency  dependent  y  and  Z^  values. 

Fig.  21  shows  the  results  for  the  complex  characteristic  impedance  for  the  1  mil 
center  conductor  copper  coax  example  for  frequencies  from  lOMHz  to  SGHz. 

The  characteristic  impedance  is  Z^  =  65.68Q -;36.9Q  at  lOMHz, 

Z^  =  58.23Q  -  y  21.0Q  at  20MHz,  Z^  =  55.0Q  -J  9.4Q  at  50MHz  and  drops  to 
Z^  =  54.13Q  -  J  5.53Q  at  lOOMHz,  still  over  8%  above  its  nominal  Z^  (nominal)  = 
50.00Q  value.  At  IGHz,  Z^  =  51.4730  -j  1.560  fairly  close  to  nominal.  For 

smaller  conductor  sizes,  as  seen  in  Fig.  22  for  the  case  of  a  1/2  mil  diameter  center 
conductor  coax  (all  dimensions  2x  smaller  than  the  1  mil  case,  so  Z^  (nominal)  still 
is  50.000),  the  effect  is  more  dramatic. 
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For  the  1/2  mil  copper  coax,  =  106.1Q -7  91.13Q  at  lOMHz, 
Z^  =  80.83Q-;59.84n  at  20MHz,  Z^  =  62.57Q -;31.0Q  at  50MHz  and 
Z^  =  57.0Q  -j  17.3Q  at  /  =  lOOMHz  (|Z^|  =  59.54Q  is  still  nearly  20%  above 
nominal  at  lOOMHz).  Even  at  IGHz,  Z^  =  52.9  -j33Q  ,  |Z^|  =  53.0Q  is  6% 

above  nominal.  Considering  that  most  analysis  methods  assume  Z^  to  be  constant, 
it  is  clear  that  problems  should  be  expected  for  these  for  small  line  dimensions. 


1.5.  Frequency  Response  of  Small  Open  -  Circuited  Coaxial 
Lines  Driven  With  a  Source  Impedance  Equal  to  their 
Nominal  Characteristic  Impedance 

As  discussed  in  the  previous  section,  we  are  up  to  the  point  in  the  frequency  domain 
calculation  where  we  have  the  values  for  the  complex  propagation  constant,  y ,  and 
complex  characteristics  impedance,  Z^,  versus  frequency.  Referring  to  the  circuit  of 
Fig.  1,  this  allows  us  to ,  given  and  Z^^  use  Eqs.  32  (or  33)  to  get  Zj  given  the 
load  impedance  Z2,  Eq.  34  to  get  V^,  and  Eq.  35  to  get  the  load  voltage,  V2.  Here  we 
will  take  the  case  of  the  Z^  (nominal)  =  50.00Q  impedance  copper  coax  lines,  with 
a  source  resistance  of  =  50.00Q  (Z,„  =  50.00D  +  j  O.OD).  We  wiU  take  the  case 

of  an  open-circuited  line  (a  CMOS-like  matched  source,  open  load  confrguration), 
Zj  =  00  ,  so  that  we  use  Eq.  33  to  get  Zj.  We  take  an  input  signal  of  =  1.0  +  j 

0.0  volts  at  all  frequencies,  in  order  that  the  magnitude  of  V2(f)  will  represent  the 
line  response  directly.  Fig.  23  shows  the  magnitude  of  KjCf)  versus  frequency  from 
lOMHz  to  7GHz  for  an  /  =  12"  length  (/^^  =  2ns)  of  1  mil  center  conductor 
diameter  copper  coax,  =  3.90,  source  terminated  with  open  load.  This  case 

corresponds  to  an  end  to  end  dc  resistance  of  about  I  =  12.2D,  which  would  be 
judged  adequate  for  digital  interconnects.  The  magnitude  of  K2  (/*)  increases  from 
l.OOOOK  at/  =  0  (dc),  to  1.0001  IK  at  20MHz,  holding  at  0.9987K  at  40MHz  before 
dropping  io\V  2if)\  -  0.95K  at  lOOMHz.  This  unexpected  rise  (from  a  nominally 
source  terminated  line,  which  should  give  a  smooth  fall)  is  part  of  an  oscillatory 
fluctuation  on  the  response  characteristic,  as  seen  in  Fig.  23.  The  first  dip  in  this 
response  is  at/  =  225MHz,  the  second  peak  at  about  300MHz,  the  second  dip  about 
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470MHz,  etc.  The  3dB  bandwidth  is  about  860MHz,  and  the  6dB  bandwidth  is 
about  3337  GHz. 

Fig.  24  shows  the  same  V2(f)  versus  /  curve  for  a  12"  length  of  1/2  mil  diameter 
center  conductor  coax.  Reducing  the  coaxial  line  conductor  dimensions  by  2x 
essentially  quadruples  the  loss  in  the  low  frequency  region,  (the  end  to  end  line 
resistance,  I  =  48.8Q  in  this  case),  and  should  double  the  loss  at  high 
frequencies.  Here,  j  Fj  (^)  |  =  0.995at  lOMHz,  0.98  at  20MHz,  0.959  at  30MHz  and 
is  0.90  at  50MHz,  falling  to  0.75  at  lOOMHz.  The  3  dB  bandwidth  is  about  120MHz 
and  the  6  dB  bandwidth  about  775MHz.  In  spite  of  the  steep  falloff  of  |  K2  (f)  |  for 
the  1/2  mil  coax  case,  the  ripple  noted  for  the  1  mil  case  in  Fig.  23  is  clearly  visible  in 
Fig.  24. 

The  reason  for  this  ripple  in  the  frequency  response  characteristics  is  the  fact  that 
the  actual  complex  line  impedances,  Z^,  as  seen  in  Figs.  21  and  22  for  these  two 
cases,  are  substantially  above  the  nominal  50.00Q  impedance  to  which  was 
matched.  For  example,  anR^  =  0,  skin  effect  only  calculation  using  Eqs.  9  and  10 
instead  of  Eqs.  40  and  41  enhances  the  ripple  ( |  K2 1  going  to  about  1.022  at  61MHz, 
for  example),  as  seen  in  Fig.  25.  Since  this  R^  =  0  case  seriously  overestimates 
at  low  frequencies,  the  corresponding  values  should  be  too  high  as  well.  On  the 
other  hand,  when  the  value  of  is  artificially  forced  to  Z^  =  50.00  +  j  O.OQ  in  the 

calculation,  instead  of  using  Eq.  28,  the  resulting  \V2(f)\  vs. /curve  (see  Fig.  26) 
falls  monotonically,  with  no  noticeable  sign  of  ripple  as  in  Figs.  23  and  24.  It  can  be 
concluded  from  this  that  the  variations  of  Z^  with /can  have  significant  effects  on 
real  MCM  interconnect  signal  quality. 
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1.6.  Application  of  Numerical  Inverse  Fourier  Transform  To 
Obtain  Time  Domain  Waveforms 

In  the  general  case  of  the  application  of  Fourier  transform  techniques  to  obtaining 
time  domain  signal  waveforms,  FFT/iFFT  pairing  is  used  to  obtain  the  convolution 
of  the  input  signal  waveform  with  the  line  response.  The  input  waveform  is 
transformed  into  the  frequency  domain  where  it  becomes  the  generator  in  Fig.  1. 
Since  the  output  of  the  transmission  line  calculation,  V2,  is  proportional  to  (and 
correspondingly  tracks  its  phase),  this  amounts  to  multiplication  in  the  frequency 
domain,  which  is  precisely  how  time  domain  convolution  is  accomplished  using  the 
FFT/iFFT  operation.  Hence,  the  iFFT  of  K2  (f  )  back  into  the  time  domain 
represents  the  convolution  of  the  input  waveform  with  the  impulse  response 
function  of  the  linear  system  (transmission  line  circuit).  While  we  may  be  utilizing 
this  software  for  this  purpose  later,  at  this  point  most  of  the  effort  has  been  on 
obtaining  the  impulse  response  of  the  lossy  transmission  line  circuit  using  the  iFFT, 
and  then  doing  a  simple  trapezoidal  rule  integration  of  the  impulse  response  result 
to  get  the  step  response. 

One  advantage  of  using  a  transformed  waveform  for  is  that  the  finite  risetime  of 
the  waveform  will  tend  to  reduce  the  amplitude  of  the  high  frequency  components 
in  V2  (f).  This  would  be  helpful,  since  if  the  transmission  line  segment  is  not  lossy 
enough,  it  may  not  be  practical  to  do  enough  points  in  the  iFFT  array  to  keep 
aliasing  of  high  frequency  components  (beyond  the  Nyguist  limit)  from  interfering 
with  the  accuracy  of  the  calculation.  Early  on,  I  successfully  used  a  step  function 
input  Oust  divide  the  K2(/')  by  j  a  to  do  this)  for  the  calculation.  While  the  l/j  (o 
factor  helps  reduce  the  high  frequency  spillover,  it  blows  up  at  zero  frequency  and 
the  error  due  to  ffnite  frequency  step  size  at  the  low  frequency  end  (between  dc  and 
the  first  ac  point),  leads  to  both  a  dc  offset  and  a  major  slope  error  which  must  be 
corrected  out  to  give  useful  results.  Hence,  I  decided  to  use  a  simple  Vj^  =  (1,0) 
frequency  domain  input  (corresponding  to  a  d  -  function  at  t  =  0  time  domain 
input),  which  makes  K2  (r)  1  the  iFFT  of  K2  (f),  the  impulse  response  of  the  system. 
If  everything  is  done  correctly,  the  iFFT  of  the  complex  K2  (f )  waveform  should 
have  a  real  part  (the  time  domain  impulse  response  waveform),  but  no  imaginary 
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part  This  was  the  case  for  the  analysis  here.  To  insure  that  the  time  step-frequency 
step  conversions  and  the  iFFT  waveform  normalization,  a  standard  Fourier 
transform  pair  (critically  damped  exponential  pulse)  was  evaluated  and  found  to 
give  very  accurate  results.  The  simple  2-point  trapezoidal  rule  integration  was 
selected  over  3-point  Simpson’s  rule  because  of  its  even  weighing  of  all  points  in  the 
integral.  Since  there  tends  to  be  some  "ping  -  ponging"  (alternate  +  and  -  errors)  in 
the  iFFT  impulse  response  due  to  feedthru  which  tends  to  average  out,  it  is 
important  that  the  integral  method  give  an  even  weighing  to  all  points.  The  results 
with  the  simple  trapezoidal  rule  integration  seem  fine,  with  step  heights  generally 
very  precisely  1.0000. 

The  case  of  the  /  =  12  inches  long  1  mil  center  conductor  coax  calculation  is  shown 
in  Figs.  27  -  29.  Fig.  27  shows  the  magnitude  of  V2(f)  versus  /  at  a  frequency  of 
61.035 15625MHz  /  point  over  the  frequency  range  off  =  0  to  25GHz.  This 
magnitude  data  is  extracted  from  the  complex  V2(f)  array  from  which  the  iFFT  is 
taken.  It  is  necessary  to  go  to  such  high  frequencies  in  order  to  insure  that  the 
magnitude  of  V2(f)  is  adequately  small  above  this  upper  frequency  limit  The 
results  of  this  8192  point  iFFT,  V2it ) ,  is  shown  as  the  impulse  response  in  Fig.  29 
(the  real  part  is  shown,  the  imaginary  part  is  zero).  This  may  be  interpreted  as  a 
voltage  slew  rate  (dVIdt),  The  peak  value  in  Fig.  29  is  a  slew  rate  of  about  20.5  volts 
/  nanosecond.  The  arrival  time  is  correct  for  an  /  =  03048  meter  line  with  =  3.90. 

Figure  29  shows  the  step  response,  V  2int  obtained  by  integrating  the  real  part  of  the 
iFFT  result,  V2(t ).  The  waveform,  with  no  correction  of  any  kind  (for  offset  or  slope 
gives  the  proper  zero  response  before  fpd  ~2ns,  a  reasonable  looking  rise,  and  an 

ultimate  value  of  l.OOOOK  It  does  have  a  slight  overshoot  "hump"  at  r  =  6ns 
(2  tpj  past  arrival),  going  to  a  peak  amplitude  of  1.021845K  at  6.084ns.  The  rise 
curve  is  03%  at  2.01ns,  2.2%  at  2.012ns,  9.8%  at  2.016ns,  20.75%  at  2.022ns,  40.13% 
at  2.040ns,  50.1%  at  2.058ns,  59.92%  at  2.086ns,  69.95%  at  2.158ns,  79.98%  at 
2332ns,  85%  at  2.53ns,  90%  at  2.910ns,  and  95%  at  3.647ns.  Past  the  "hump",  the 
step  amplitude  drops  to  0,999940298K  at  13.322ns,  ending  at  0.99999998 IK  at 
16.328ns  in  the  8192  point  array  (2.0  ps/step).  The  "hump"  response  is  real,  the 
result  of  the  frequency  dependent  Z^,  corresponding  to  the  "ripple"  seen  in  V2(f ).  In 
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the  case  of  the  =  0,  skin  effect  only  analysis  shown  in  Figs.  30  and  31  for  the 
same  /  =  12"  1  mil  center  conductor  diameter  coax,  the  result  is  a  more  pronounced 
"hump"  with  nearly  4%  overshoot,  corresponding  to  the  higher  ripple  in  V2(f)  seen 
in  Fig.  25  for  this  =  0  case.  On  the  other  hand,  as  shown  in  Fig.  32  for  the  case 
where  was  forced  to  be  Z^  =  SO.OOQ  +  j  0.0O(as  shown  in  Fig.  26),  there  is  no 

hump  and  a  remarkably  erfc-like  behavior  (as  suggested  by  Eqs.  14-18)  is  observed. 
This  would  appear  to  point  to  the  "hump"  as  being  a  real  consequence  of  the 
frequency  dependence  of  Z^  in  these  small  lines,  and  not  some  calculational  artifact. 

Figures  33-35  show  the  case  for  the  1/2  mil  diameter  center  conductor  coax,  also  /  = 
12  inches  long.  Because  of  the  higher  alternation  in  the  smaller  center  conductor, 
little  signal  output  is  seen  beyond  lOGHz  in  Fig.  33.  Correspondingly,  the  maximum 
voltage  slew  rate  in  Fig.  34  is  smaller,  less  than  5  volts/ns.  The  step  response  for  this 
I  =  48.8Q  case  shows  the  same  kind  of  pronounced  R^  sloping  segment  effects 

discussed  at  length  in  last  month’s  report,  and  seen,  for  example  in  Fig.  12. 

In  summary,  the  results  of  this  full-blown  frequency  domain  /  iFFT  intercoimect  line 
simulation  appears  to  be  physically  meaningful  and  ofrer  important  insights  into  line 
behavior. 
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Rac  and  Lint  vs.  frequency  Calculated  using  Kelvin  Functions  for  1  mil  diameter  50Q  Copper  Coax 
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Section  2.0 

Commercialization  of  High  Temperature 
Superconducting  Technology 


2.1.  Background 

For  the  quarter  ending  June  1992  much  of  the  emphasis  was  directed  at  working  out 
a  cryoelectronic  plan  that  could  be  sold  to  DARPA  management  and  one  that  would 
culminate  in  an  industry  driven  by  this  technology.  It  is  now  accepted  that  while  the 
high  end  i^stem  market  may  be  the  first  to  start  designing  into  ciyoelectronic,  the 
stable  part  of  the  industiy  has  got  to  come  from  the  high  part  of  the  low  end  such  as 
the  workstation.  That  is  where  the  volume  market  will  be  and  that  is  where  we  think 
we  can  make  a  significant  impact 

When  we  began  this  project  we  were  thinking  of  it  as  a  long  term  program  waiting 
for  breakthroughs  in  HTSC  materials  and  MCM  processes.  This  approach  has  now 
been  expanded  significantly.  We  have  recommended  to  DARPA  that  HTSC-MCM 
should  be  considered  as  part  of  a  bigger  ciyoelectronic  technology  which  goes 
b^ond  TT’K  operation  and  start  thinking  about  a  cascaded  approach  which  begins 
at  around  20(f  K  and  works  its  way  down  to  '/TK.  In  view  of  the  fact  that  future  IC 
foundries  will  be  very  costly,  in  excess  of  $1B  per  foundry,  and  in  view  of  the  fact  that 
larger  IC’s  may  not  have  the  performance  that  has  been  anticipated,  it  may  be  that 
low  temperature  IC  operations  could  bring  in  a  much  more  favorable 
price-performance  characteristic  in  a  very  near  short  time  using  a  modified  and 
optimized  IC  process  directed  for  low  temperature  operations.  To  our  pleasant 
surprise  we  find  the  same  type  of  thinking  in  industiy  be  it  high  end  or  low  end 
computer  systems. 
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Durinfj  the  last  few  months  a  number  of  meetings  were  held  at  SUN,  STl,  Carrier, 
Cray,  £-Systems,  Datamax  and  DEC.  All  of  these  companies  now  seem  to  be 
converging  towards  a  cryoelectronic  technology  thrust  that  can  meet  oin  business 
criteria  close  to  an  affordable  budget. 

The  bottom  line  is  this:  a  merchant  cryoelectronic  mfrastructure  led  by  a  company 
such  as  Sn  needs  to  be  developed  to  service  the  ^sterns  industry.  This  company 
would  want  to  supply  ciyo-subassembiles  to  system  houses.  They  would  work  IC 
foundries  to  optimize  IC  for  cryoelectronics  and  set  up  the  mechanics  to  market 
them  to  the  industry.  They  would  work  with  the  CAD  tools  industry  to  ensure  a 
complete  suite  for  the  cryoelectronic  industry  is  available  and  that  equipment  for 
packaging  and  testing  the  ICs  can  be  developed.  Most  importantly  they  would  work 
to  see  that  a  cost  effective  highly  reliable  cryogenics  be  available  within  two  years. 

Extensive  discussions  of  this  approach  took  place  with  SUN,  Cray,  and  DEC. 
Programs  with  these  companies  are  proprietary  and  cannot  be  discussed  in  this 
report  However,  based  on  these  discussions,  a  four  year  cryoelectronics  plan  was 
developed  and  submitted  to  E  Patten  at  DARPA.  In  this  plan,  details  were  given 
about  the  requirements  for  infrastructure  technology  elements  development, 
applications  both  commercial  and  military,  required  budgets,  milestones,  and 
timelines.  This  report  can  be  accessed  through  Dr.  Patten. 
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